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ABSTRACT 

We present new results on the physical nature of infrared-luminous sources at 0.5 < z < 2.8 
as revealed by HST/NICMOS imaging and IRS mid-infrared spectroscopy. Our sample consists of 
134 galaxies selected at 24/im with a flux of S(2Ajim) > 0.9 mJy. We find many (~ 60%) of our 
sources to possess an important bulge and/or central point source component, most of which reveal 
additional underlying structures after subtraction of a best-fit sersic (or sersic+PSF) profile. Based 
on visual inspection of the NIC2 images and their residuals, we estimate that ~ 80% of all our 
sources are mergers. We calculate lower and upper limits on the merger fraction to be 62% and 91% 
respectively. At z < 1.5, we observe objects in early (pre-coalescence) merging stages to be mostly 
disk and star formation dominated, while we find mergers to be mainly bulge-dominated and AGN- 
starburst composites during coalescence and then AGN-dominated in late stages. This is analogous 
to what is observed in local ULIRGs. At z > 1.5, we find a dramatic rise in the number of objects 
in pre-coalescence phases of merging, despite an increase in the preponderance of AGN signatures 
in their mid-IR spectra and luminosities above 10 125 L©. We further find the majority of mergers 
at those redshifts to retain a disk-dominated profile during coalescence. We conclude that, albeit 
still driven by mergers, these high-z ULIRGs are substantially different in nature from their local 
counterparts and speculate that this is likely due to their higher gas content. Finally, we observe 
obscured (rg.7 Mm > 3.36) quasars to live in faint and compact hosts and show that these are likely 
high-redshift analogs of local dense-core mergers. We find late-stage mergers to show predominantly 
unobscured AGN spectra, but do not observe other morphological classes to occupy any one specific 
region in the Tg^fim vs - PAH equivalent width (or Spoon) diagram. This suggests a high degree 
of variation in the PAH emission and silicate absorption properties of these mergers, and possibly 
throughout the merging process itself. 

Subject headings: galaxies: evolution — galaxies: high-redshift — infrared: galaxies 



1. INTRODUCTION 

Since their discovery in the IRAS all-sky survey more 
than 25 years ago, ultra-luminous infrared galaxies 
(ULIRGs, L IR > 10 12 L©) have been thought to rep- 
resent a key evolutionary link between normal galaxies 
and quasars. Early evidence suggested that they were 
young quasars fed and obscured by large amounts of gas 
and dust funneled towards the center of a mer ger rem- 
nant of two gas-rich spirals ( Sanders et al.|1 988). Intense 
star formation rather than black hole accretion was later 
shown to be the primary energy source of mo st ULIRGs 
flRigopoulou et al.||1996| |Genzel et aL]|1998| ), but black 
hole accretion has remained the dominant mechanism at 



higher luminosities (L IR > lO^ 2 - 3 !^; |Lutz et al.| 1998; 
Veill eux et aTJ[l999l |Tran et al.|[200l| |Farrah et al.| 2003 ; 



Veilleux et al. 2009b| ). Ground-based as well as high 
resolution HS'l imaging, meanwhile, revealed that more 
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than 95% of ULIRGs originate in a merger, but that the 
ULIR G phase can also appear much before fina l coales- 
cence flMurphy et aL|l996[ [Veilleux et al.||2002[ ). 

Simulations have yielded strong support to the pic- 
ture of quasars originating in the merger of two gas- 
rich spirals after a phase of intense star formation and 
rapid black hole growth giving ri s e to t he UL IRG phe- 
nomenon. iBarnes fc Hernquist (1991) and Mihos & 



|Hernquis t (1994) showed that the merger of two equal- 
mass disk galaxies can rapidly dissipate angular momen- 
tum, causing the gas to fall to the center of the galaxy 
and create a starburst of ULIRG proportions. They also 
demonstrated how the stellar component formed tails 
and streams much like the ones observed in ULIRGs. 
Refinements in hydrodynamical s imulations and the in- 
troduc tion of AGN feedback by Spri ngel fc Hernquist 
(2005) then showed how the AGN, once triggered, can 
expel remaining gas, quench star formation and become 
a true quasar w ithin a remnant elliptical galaxy ( |Di Mat-] 
|teo et al.|2005l|Hopkins et al.|2008[ ). Finally, the addition 
ol radiative transfer confirmed the exceptional infrared 
luminosity associated with the w hole event, and in p ar 



ticular with the final coalescenc e (|Jonsson e t al. 20 061 |Li 
et al.|2008[ | Younger et al |2009| [Narayanan et al.|2UlUa| T 



Although ULIRGs and quasars are extremely rare lo- 
cally (|Soifer et al.||1987 ), observations at sub- millimeter, 
mid-IR, optical and X-ray wavelengths have demon- 
strated that their number and luminosity densities in- 
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crease rapidly wit h redshift ([Chapman et al. 2QQ5[ Le 



Floc'h et al.||2005| [Richards et al.||2UU6| |Hasinger et al. 



2QQ5p . In particular, the advent of the Spitzer Space Tele- 



scope has enabled sensitive and fast imaging at 24 /im 
with MIPS, yielding the detection of a large number of 
infrared- lum inous galaxies at z < 3 (e.g. |Perez-Gonzalez 



et al.||2005| ). With such a rise in prominence, it becomes 



important for our understanding of galaxy /quasar evo 
lution to ask whether, or how many of, these numerous 
high-z ULIRGs are triggered through the same physi- 
cal mechanisms as their low-redshift counterparts, and 
whether they also represent a tr ansition towards quasars . 
Given the higher gas f ractions ([No terdaeme et al. 2009), 
star formations rates ( [Hopkins "ST Beacom 2006) as well 
as specific star formation rates ([Zamojski et al. 2007 



selection criteria. In an effort to address this question 
more fully, we have obtained HST/NICMOS imaging of a 
sample of 135 bright high-redshift 24/im-selected galax- 
ies that combines the previousl y published data for 33 
sources from our first program (|Dasyra et al.||2008|) to 
new data for 102 sources from our second progra m. 

In t his paper, we expand the analysis of Das yra et al. 
( 2008 ) and perform a systematic search of merger signa- 
tures by uncovering underlying structures, resulting in 
the discovery that ~ 80% of our objects are fueled by 
either an ongoing or recent (< 0.5 Gyr) merger event. 
We also artificially redshift local ULIRGs to quantify 
the detectability of merging signatures at high redshift 
and to create a comparison sample for our sources. We 
examine similarities and differences with local ULIRGs. 
Brinchmann & Ellis 2000) of the overall galaxy popu- We then combine our morphological results with our IRS 



lation at these redshift s, quiescent star formation is ex- 
pected to contribute i ncreasingly more to t he infrared lu- 
minosity of galaxies ( Hopkins et al.pOlO ). Certainly at 
z ~ 2, some ULIRGs nave been found toexhibit disk-like 
kinematics (|Forster Schreiber et al. 2009 Carilli et al. 



20101 iBothwell et al.||2010|). The answer as to what ex- 



tent mergers are still necessary to explain the origin of 
the infrared-luminous population at higher redshifts is, 
therefore, unclear. 

Meanwhile, ISO studies of local ULIRGs have demon- 
strated that mid-IR spectroscopy, through the resolu- 
tion of PAH emission complexes and measurement of 
their strength relative to the underlying continuum, is 
the most effective single tool for identifying which of star 
formation or AGN activity is responsible for the observed 
mid-IR radiation of an object. The launch of the Spitzer 
Space Telescope with its infrared spectrograph (IRS) sen- 
sitive to fluxes of 5 , 24 /i m ^ 1 niJy has, thus, brought a 
flurry of mid-IR spectroscopic surveys of bright, high- 
redshift 24/im galaxies ai med at addressing the origin of 
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their infrared luminosity (|Houck e t al. 2005 Yan et al. 
Weedman et al.|2006||Yan et al |2007| |b>ajma et al. 
Farrah et al.||2008[ |Dasyra et al.||2009[lDesai et al. 



2009). The results of these efforts have demonstrated 
that, unlike sub-millimeter galaxies that are primarily 
powered by star forma tion ( |Pope et al. [2008 1 iMenendez- 



Delmestre et al.| [2009), 24/im-selected objects appear to 



be more analogous to local ULIRGs in that they display 
both types of spectra (as well as various combinations 
thereof). Also in analogy to local ULIRGs, the AGN 
contribution to their mid-IR flux increases with total IR- 



lumin osity (Sajina et al. 2007 Dey et al. 2008 Desai et al. 



2009]). 

Our group carried out two such mid-IR spectroscopic 
programs of 24/im-bright galaxies, the firs t of which 
yielded spectra for 52 objec ts at z > 1 - 3 ( |Yan et aL 
|2007[|Sajina et al.|2007[|2008[ ), while the second targeted 
150 objects spanni ng a red shift range of z~ 0.3 -2.5 and 
peaking at z = 1 (Dasyra et al. 2009). Both were con- 
ducted in the Spitzer Extragalactic First Look Survey 
(XFLS). Analysis of this combined data has revealed the 
presence of an obscured AGN in > 75% of our objects. 

The answer as to whether the infrared luminosity of 
these 24/im-bright galaxies is due to mergers, however, 
has, thus far, not been conclusive. I nitial morphological 
studies have yielded mixed results ([Dasyra et al. 2008 



|Bussmann et al.| |2009 ; M elbourne et al.||2009| ), but part 
ol this variation could be due to small sample sizes and 



spectra ( |Sajina et al. 2007| |Dasyra et al." [20091 ) and the 
SED analysis ol Sajina et al.l (|2008|) and Sajina et al. 



(in preparation) to investigate tlie~link between mor- 
phology/merging and, both, the relative strength of the 
AGN and starburst components at infrared wavelengths, 
as well as the degree of obscuration. We discuss our re- 
sults in the context of other LIRG and ULIRG samples 
at both high and low redshift. Finally, we confront our 
observations with simulations, and discuss their conse- 
quences for our understanding of galaxy evolution. 

The paper is organized such that we first describe our 
sample and analysis in §2, which we follow by a discussion 
of our simulated observations of redshifted local galaxies 
in §3. We present our results in §4 and compare them to 
that of other ULIRG samples in §5. We end section 5 by 
discussing implications for the theory of galaxy evolution, 
and conclude with a summary (section 6). We use a 
ACDM cosmology with H = 70 km s" 1 Mpc -1 , Q m = 
0.3 and Q A = 0.7. 

2. OBSERVATIONS & ANALYSIS 
2.1. Sample 

Our full sample consists of 134 galaxies with 24/im 
fluxes above 0.9 mJy selected from the Extragalactic 
First Look Survey. It draws from two sepa rate pro- 
grams whose IRS sp ectra wer e presented in [Yan et aL 
Q2007| ); |Sajina et al.| ( [2007| ) and |Dasyra et al.| ( [2009| ^We 
have obtained H5T/NTC2 F16l)W (^-band) iSages for 
33 sources from our first program and 102 sources from 
our second program. We rejected one source that we sus- 
pect of being gravitationally lensed. All data were taken 
in MULTIACCUM mode, and we obtained, for each ob- 
ject, four to eight dithered images with exposure times 
ranging from 640 to 672 seconds for a total exposure time 
of 2560 to 5376 seconds per o bject. The data for t he first 
33 sources were presented in IDasyra et al.| |2008[ ) . 

In parallel with our NICMOS observations, we ac- 
quired WFPC2/F814W images of patches of the FLS 
field visible to that instrument at the time of our primary 
observations. Eight of our objects fall in one or another 
of these pointings. An additional twelve of our objects 
fall in the central 0.12 deg 2 of the FLS field that has been 
image d by HST/ACS in the F814W filter ( [Bridge et af 
2007j ). 

Galaxies in both our programs were chosen randomly 
from the entire set of 5(24/ira) > 0.9 mJy sources in the 
FLS main field using, in the first case, a color selection of 
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vf v (2^m)lvf v ^m) > 3.16 and vf v (2^m)/vf v (R) > 
10, and in the second case, a sampling rate varying ac- 
cording to their i?-band magnitude, going from 1 in 10 
for objects with 20 < R < 22 to 1 in 3 for objects with 
22 < R < 24.5, while avoiding objects already targeted 
in our first program. The strong weight put on i?-band 
faint objects ensures that a large fraction of them lie 
at high redshifts (z ~ 1-2), thus preventing our sample 
from being populated at ~ 70% by low redshift (z < 0.3) 
sources as is the case for the general S'(24/im) > 0.9 mJy 
population. Our initial cut on the 24/im/8/im color was 
made to isolate objects with eith er strong PAH emis- 
sion or a steep mid-IR continuum ( |Yan et al. 2007), but 
this bias is compensated for in our combined sample as 
shown in Figure [I] Figure [I] illustrates the distribution of 
our sample in S 2 4 Mm /^ M m vs. S 2 ^ m / 'So. 64 /im color-color 
space in comparison to the general population of bright 
24/im-sources in the FLS field, and Figure [2] shows its 
distribution in redshift. 
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Fig. 1. — Color-color plot comparing objects from our sample (red 
and magenta squares) to the general population of bright 24/im- 
sources (dots) in the FLS field. Magenta represents objects drawn 
from our first program and red from our second. Dashed and solid 
lines represent selection criteria used in | Dasyra et al.| ( |2008| and 
|Bussmann et al.] ( |2009| > respectively. 



In this paper, we utilize redshift measurements ex- 
tracted from t he IR S data using the method described in 
Sajina et al.| ( |2QQ7 ). When available, we use more pre- 



cise redshifts derived from optical or near- infr ared spec- 
tra instead (~ 20 objects; Sajina et al.||2QQ8| ). We miss 
redshift estimates for twelve of our galaxies. These are 
objects that show absolutely featureless mid-IR spectra, 
and that have not been observed spectroscopically in the 
optical/NIR either. Among the ones we do have lines 
and features for, our lowest redshift object is found to be 
at z = 0.24 and our highest one is at z = 3.48, but most 
of our sources fall in the range 0.5-2.8 with a main peak 
at z ~ 1, and a secondary peak at z ~ 2. 

Using ancill ary far-infrared, su b-millimeter and radio 
observations, |Sajina et ah (2008) and Sajina et al. (in 
preparation) were able to constrain the SED and de- 
rive full (3-1000/im) infrared luminosities for 112 ob- 
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Fig. 2. — Reshift distribution of our sample. We miss redshift 
estimates for twelve of our objects. 



jects in our sample. Their analysis tells us that most 
of our sources possess luminosities ranging from 10 11,5 to 
lO 13 I/ , and that, because we have a flux-limited sample, 
higher redshift objects also possess higher luminosities, 
such that most objects at z > 1 have an Lir > 10 12 L© 
and nearly all objects at z > 1.5 have luminosities above 
10 12,5 L©. We discuss the relation betw een redshift and 
luminosity in more detail in section |4~4| 



2.2. Control Sample 

We use galaxies present around our 24/im-selected 
sources in the NICMOS images as a control sample for 
our morphological classification. Figure [3] shows the 
number of galaxies we find in our control sample in bins 
of 0.4 magnitude. For comparison, we also plot the num- 
ber of galaxies we would expect to find in our ens emble 
of images b ased solely on the nu mber counts of [Chen 
et al.|fl2002l ) an d |Yan et afl (|l 998). Most r ecent studies 
such as those of [Metcalfe et al. ( 2006]) and |Retzlaff et ah] 
(2010) arrive at nearly identical numbers. Our control 
sample shows an excess at H ~ 18-20, where our sam- 
ple of 24/im-selected galaxies peaks. We interpret this as 
representative of an over-density in the environment of 
24/im-sources. The remainder of our galaxy counts, on 
the other hand, agrees remarkably well with published 
values. 

It is important to note that in counting galaxies, we 
exclude all parts associated with the 24/im source, and 
retain only the random component (that is the surround- 
ings). There is ambiguity, in a handful of cases, as to 
whether some distant objects are associated with their 
neighboring 24/im-galaxy or not. We chose to include 
only half of them in our number counts, but since there is 
no more than two such cases in any given bin (many have 
none), this never amounts to more than one count. All 
are folded into our uncertainties though. In our account- 
ing, we also exclude all point sources with magnitudes 
brighter than H = 22 in order to avoid contamination 
by stars. Such a cut, however, also rejects extra-galactic 
point sources (QSOs). It is therefore likely that we are 
slightly underestimating our counts, though again, this 



4 



Zamojski et al. 



100.0 



Control somple 
Chen et ol. (2002) 
Yon et ol. (1998) 
Fit 

24/zm somple 



6 



10.0 




.o 
E 



H mag 

Fig. 3. — Number of galaxies per bin of 0.4 magnitude in our con- 
trol sample (blue curve and points). For compa rison, the expected 
nu mber of galaxies f rom the number counts of |Chen et al.| (2002) 
and |Yan et al.| |l998 ) are shown in open circles and triangles respec- 
tively, along with a fit to those points (red curve). Uncertainties 
are shown only when larger than the symbol size. The histogram 
shows the magnitude distribution of our 24/ira-selected galaxies. 
We observe an excess in our number counts at magnitudes where 
our 24/ira-selected sample peaks. 



can only be a small effect. 

2.3. Data Reduction 

We used a modified version of the standard NICMOS 
reduction pipeline created by V. Fadeyev, and that prop- 
erly corrects for the effect of cosmic ray s on neighbor- 
i ng p ixels. This routine is described in Fade yev et al.| 
(2006). Following the standard procedures, we used the 
cainica, nicpipe, and biaseq routines to perform the ba- 
sic reduction steps, with bias equalization, on each of 
our images. We then ran pedsub to remove the "pedestal 
effect "between the four quadrant of the images. For im- 
ages affected by the South Atlantic Anomaly, we then ran 
saaclean to remove the cosmic rays imprint that persists 
for a certain time after the telescope has left the SAA- 
region. We then ran the pedsky routine to perform a first 
sky subtraction. 

These first standard steps correct most instrumental 
signatures and bring the background to an average of 
zero. However, at the end of this process remaining resid- 
ual patterns are often present in many images. Our next 
step was, thus, to create a "supersky "image by median 
averaging images that still carried a strong residual pat- 
tern in their background, and smoothing that median im- 
age with a ring median filter with inner and outer radii 
of two and ten pixels respectively. At the end, we used 
25 images to construct this supersky image. We then, for 
each image, found the best-fit factor by which to multiply 
the supersky such that, when subtracted from the image, 
minimizes the dispersion in that image; and proceeded 
to subtract the scaled supersky. Finally, we performed a 
background column subtraction on images that presented 
vertical features usually caused by electronic ringing and 
streaking of bright sources. 

We then combined all the dithered images of each ob- 



ject into a final mosaic using the dither package part 
of the STSDAS external package in IRAF. We first driz- 
zled each final NICMOS image onto separate, but aligned 
grids, sampled at half the size of a NIC2 pixel, using a 
square kernel (or "drop ") of 0.8 times the NIC2 pixel 
size. We then combined the drizzled images into one 
final image using a 3<r rejection around the median. 

The surface brightness detection limit in our final im- 
ages is ~ 23.9 mag arcsec -2 , given our choice of ex- 
traction parameters, which consists of smoothing with 
a 3 x 3 pixels pyramidal kernel and requiring that at 
least 36 pixels above 1.5<j of the noise, whose RMS is 
22.5 mag arcsec -2 , be connected to one another. This 
detection limit corresponds well to the surface bright- 
ness we are able to pick up by eye. Accounting for sur- 
face brightness dimming, it is equivalent to an intrinsic 
surface brightness of \i — 20.9 mag arcsec -2 at redshift 
one and \i = 19.6 mag arcsec -2 at z = 1.7. 

2.4. Profile Fitting 



We used the GALFIT package flPeng et al. |2002[ ) to 
perform a series of fits to the profile ol all of our ob- 
jects. We first fitted our objects with a one-component 
sersic profile with free sersic index. For those objects 
with a sersic index, n, greater than 1 that were suf- 
ficiently extended and whose residuals (or image) sug- 
gested the presence of a disk, we repeated the fit with a 
two-component bulge (de Vaucouleurs) plus disk (expo- 
nential) profile. We use this bulge-to-disk decomposition 
to classify our objects into disk-dominated, B/D < 1, 
and bulge-dominated, B/D > 1, galaxies. 

Because many of our objects consist of a bulge sur- 
rounded by an envelope of tidal streams, the addition 
of a disk component to the fit often results in a large 
fraction of the flux of those objects to be assigned to 
that component. This is unphysical since these tidal en- 
velopes do not follow an exponential profile. In those 
cases, we instead assume that the bulge is well repre- 
sented by the sersic profile and that the disk component 
is negligible. Objects that are too small to be meaning- 
fully decomposed into bulge and disk components, on the 
other hand, we classified using the value of their sersic in- 
dex only. Comparison of the fitted sersic index with the 
bulge-to-disk ratio of the more extended objects in our 
sample indicates that a5/D of 1 roughly corresponds to 
a sersic index of 2.3. We therefore adopt n = 2.3 as the 
limit between bulge and disk-dominated galaxies. 

Many (~ 20%) of our objects display a bright PSF pro- 
file at their centers, suggestive of the presence of a sig- 
nificant point source. Most of them turned out very high 
sersic indices in their one-component fits. We, therefore, 
re- fit those objects for which the one-component, sersic 
profile, yielded an n > 4.5, with a two-component, PSF 
plus sersic, profile. Four galaxies that looked like they 
might have a point source turned out indices of n ~ 4, 
only with small effective radii, and were therefore not 
re- fit. Some of our largest galaxies, we were able to suc- 
cessfully fit with a three-component, PSF+bulge+disk, 
model. We also added a PSF component to the fit of 
some of our galaxies when the one-component and two- 
component fits showed significant residuals and those 
residuals largely disappeared when the PSF component 
was added. In all cases, we always constrained the mag- 
nitude of the PSF component (when present) using only 
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a PSF+sersic fit, and then proceeded to fit the bulge 
and disk components, when possible, maintaining the 
PSF component fixed. We then used the bulge-to-disk 
decomposition or sersic index of the host to determine, 
when possible, whether these objects were disk or bulge- 
dominated. For a number of objects, however, the host 
galaxy was too faint compared to the central PSF and too 
small for the sersic profile to be reasonably constrained. 
The profile of those objects was therefore left ambiguous. 

Many of our objects have companions nearby. In such 
cases we always fitted the main optical counterpart to our 
24/im source and its companions simultaneously. The 
classification of all our objects as either bulge or disk- 
dominated can be found in Appendix B, along with their 
magnitude, size, and morphological class. Morphological 
classes are described in the following section. 

2.5. Morphological Classification 

Locally, ~ 99% of ULIRGs visuall y appear to be on - 
going mergers or merger remnants ( Veilleux et al.|2 QQ2). 
Analogously, a large number of our objects show signs 
(or possible signs) of merging. In order to establish a 
solid comparison with local U LIRGs we chose to u se the 
same classification scheme as Veill eux et al. J 2QQ2 ). This 
scheme, originally proposed by |Surace ( |l99oj) , classifies 
objects according to their stage ol merging, following a 
merging sequence. It is useful in that it allows to study 
changes in galaxy properties as a function of merging 
phase, and thus uncover how these properties change, 
statistically, througho ut a merger event. We present 
this schem e in section 2.5.1 We then introduce, in sec- 
tion (2T5T2I confidence classes which we use to group ob- 



jects according to how confident we are about their merg- 
ing nature. We show examples of objects from each of 
our confidence classes. Finally, we briefly ment ion auto- 
mated classification techniques in section [2. 5. 3| 

2.5.1. Classification Scheme 



In th is p aper, we use the classi fication scheme of |Surace] 
(19981 and lVeilleux et gil.] ( |2QQ2| ) which partitions objects 
into the following classes , following the progression of a 
merger: 

I. First Approach or Pair - This category refers to 
galaxy pairs that have not passed through each other 
yet. In this phase the two objects still retain their own 
characteristic morphology (albeit maybe with some per- 
turbations). 

II. First contact - At this stage, the two objects over- 
lap, but tidal streams and debris have not formed yet. 

III. Pre-merger - After their first encounter, the two 
galaxies will usually emerge again (except in very par- 
ticular low-speed, co-aligned configurations where they 
could merge during first pass), but this time tidal tails 
and bridges will have formed, and their respective mor- 
phology will be highly disrupted. Double-nucleus sys- 
tems showing tails, bridges or morphologies suggestive 
of a recent encounter fall in this category. 

IV. Advanced merger- At this stage, the two nuclei are 
undistinguishable meaning they have either coalesced or 
are on the verge of doing so. Trails of stars that have 
been tidally stripped or other debris from the merging 
process are readily visible after, and sometimes even be- 
fore, subtraction of the main, smooth component of the 



galaxy. This is also a phase of rapid bulge growth. At 
high redshifts, tidal features can sometimes appear as 
extended lopsided disks. In those cases, we require that 
the main component be a bulge in order for the object 
to be classified as an advanced merger, so that not to be 
contaminated by disturbed, asymmetr ic or clumpy disks 
whose ori gin might not lie in a merger ( Forster Schreiber 
eTaLl|201Q ). 



V. Old merger - In this phase, tidal tails and streams 
have faded out, but the galaxy still shows signs of 
past events in the form of residual asymmetries and/or 
clumps, that are detectable through subtraction of the 
smooth profile. Because signatures typically associated 
with mergers have disappeared from the objects classified 
in this category, they cannot be known for sure to have 
experienced a recent merger. We require them, however, 
to be bulge or PSF-dominated, so that their observed un- 
derlying structure, which must be important since still 
detected at z > 1, when put together with a bulge-like 
galaxy, is most easily explained through the merger sce- 
nario. 



I Veilleux et al. ( 2QQ2| further split phase III into close 
and wide binaries! We do not make that distinction, 
although we would have about equal numbers in each 
class. The largest separation we o bserve among our pre- 
mergers is 35 kpc (in MIPS 298). [Veilleux et aT] ( |2QQ2 ) 
also split phase IV into diffuse and compact mergers, but 
the lower physical resolution of our data does not allow 
us to make that distinction either. On the other hand, 
we find four triplets in our sample and have a category 
for them. However, because stages II, III and triplets all 
represent objects that are in the process of merging, but 
have not yet coalesced, we usually refer to them collec- 
tively as early mergers. 
We note tha t , beca use t he imagery of local U LIRGs 



by |Kim et aL 
much higher p 



(|2002j) and |Veilleux et al.| ( 2QQ2[ ) has a 
ihysical resolution and much higher surface 
brightness sensitivity than we have for our high-redshift 
sample, their phase IV covers a larger part of the merging 
process than our phase IV, as we lose the streams and 
tails that characterize that phase earlier than they do. 
Furthermore, it is impossible for us to detect the resid- 
ual artifacts of a merger to the level they do to define 
their phase V. In our data, many of those objects would 
look like regular bulges (see section 3). Our phase V ob- 
jects, assuming they are truly mergers, should rather be 
thought of as a mix of late phase IV and early phase V 
objects in local terms, a kind of stage 4.5. 
Although many of our objects are mergers (cf. 



4.1) 



our sample does also contain a number of regular-looking 
unperturbed galaxies. On top of the six (including 
triplets) merger categories described above, we therefore 
also have: face- on and edge- on spirals when either spiral 
arms or an edge-on disk are readily distinguishable, and 
regular bulges when no particular features are visible and 
the galaxy has a bulge-to-disk ratio B/D > 1, or a sersic 
index n > 2.3. We also have one object that appears as 
a pure point source in our data. 

Lastly, we have a number of objects in our sample that 
we label faint & compact, because they are best fit with 
a disk-ish (n ~ 1-2) profile, yet are much smaller and 
fainter than regular spiraqj and show no signs of either 
spiral arms or of an extended disk. Most of them show, 
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instead, stage V type of residuals. We show in section 3 
that some local ULIRGs can indeed appear as faint & 
compact at high redshift. They are typically phase IV 
objects that possess diffuse tidal tails too faint to be vis- 
ible at high redshift s, and dense cores often intersected 
by dust lanes, rather than fully formed bulges. Because 
of that fact, we often group these faint & compact ob- 
jects with phase IV galaxies proper and refer to them as 
coalescence objects. 

2.5.2. Confidence classes 

Although there is inherent uncertainty associated with 
morphological classification, some objects can be classi- 
fied with more confidence than others and, consequently, 
some results are more secure than others. It is, thus, 
useful to distinguish these degrees of certainty. We do so 
by assigning to all of our objects a confidence level with 
which their merging nature can be inferred. Although 
these levels of confidence can themselves be somewhat 
subjective, we find them, from a second, independent re- 
classification by another co-author (L. Yan), to vary by 
no more than one confidence level in 80% of cases, show- 
ing that they are, in practice, reasonably well-defined. 
We look below at what we can learn from each of these 
confidence classes. 

We call the most secure objects, those that can be im- 
mediately and unarguably identified as mergers at first 
sight, category 1. These objects consist exclusively of 
two overlapping or connected galaxies that show tidal 
tails and streams of their interaction, as these are the 
only kind whose merging nature we find can be estab- 
lished without a doubt from visual examination alone. 
A local analog would be The Antennae galaxy. Unfortu- 
nately, we find only five objects that fit that description, 
in our sample. They are all shown in Figure |1J Besides 
containing a low number of galaxies, this first category 
further includes only objects in stage III of the merging 
process or triplets. It thus becomes immediately clear 
that in order to make any progress, we need to move 
beyond this realm of absolute certainty. 

We thus turn to the next best objects which we call 
highly probable mergers. We include in this second cate- 
gory objects whose morphology is strongly suggestive of 
a merger event, but cannot entirely exclude other pos- 
sibilities^] An example of which would be close pairs, 
since without kinematics, their merging nature cannot 
be asserted with 100%, but the proximity of the galaxies 
to one another (typically < 20 kpc in projection) and 
their comparable size strongly suggests that they are in 
the process of merging. Another example would be ad- 
vanced (singly nucleated) mergers with highly distorted 
morphologies and strong tidal tails, as one could imagine 
features like that arising in a process such as a high-speed 
encounter, even though they are far more typical of merg- 
ers. Figure [5] shows examples of different highly probable 
mergers. 

7 Faint Sz compact objects have typical half-light radii and ap- 
parent magnitudes of 7*1/2 ~ 0.15-0.35" and mn > 19.6, whereas 
spiral galaxies in our sample have values of 7*1/2 ~ 0.4-0.6" and 
17.7 < m H < 19.0. 

9 Because the exact probabilities for an object in each of our con- 
fidence classes to be a merger are unknown, we are forced to rely on 
approximate terminology such as "highly probable" , "strongly" and 
"far more" . 



By including highly probable mergers, we more than 
quintuple our number of mergers to a total of 27, and 
add very little uncertainty in the process. We also span 
all of the major stages of merging (four close pairs [stage 
I], one phase II object, seven more stage III mergers, 
one more triplet, and nine advanced mergers [phase IV]). 
Categories 1 and 2 together, however, still represent only 
20% of our sample. 

Categories 3 and 4, illustrated in Figures [6] and [7] re- 
spectively, include, on the other hand, 35 and 33 ob- 
jects each. These objects also show merger signatures, 
although these progressively go down in the strength 
and/or shape of tidal features, the distance and/or con- 
trast ratio between the primary object and the compan- 
ion, and in the overall connectivity of the system. The 
merger origin of these objects is therefore increasingly 
more uncertain. It is nevertheless clear that they must 
contribute an important number of mergers, although we 



postpone more specific estimates until section 4.1 



We then also have category 5 objects. These are ob- 
jects that no longer show companions or tidal features 
characteristic of mergers, but do show residuals and 
asymmetries after subtraction of the smooth component. 
The majority of our phase V objects (the late mergers) 
and of our faint & compact objects fall in this category, 
examples of which are shown in Figure [§J Some local 
ULIRGs, that are clear merger remnants, also show sim- 
ilar features when redshifted to z > 1, as demonstrated 
in section 3. Therefore, although asymmetries can arise 
in numerous situations, we must still consider it proba- 
ble that these objects have their origin in a merger event. 
This is true for phase V objects as their typically high- 
sersic profiles combined with residuals bright enough to 
be detected at z ~ 1-2 offers few alternatives, but also 
of Faint & Compact objects as argued in § 3. 

Finally, category 6 consists of elliptical galaxies that 
are well-fit with a smooth sersic profile and show no 
residuals. The reason why one might choose to count 
some of these objects towards the merger fraction orig- 
inates from the fact that, as demonstrated in section 3, 
late mergers often lose their merging signatures in the 
redshifting process. Regular ellipticals could, therefore, 
be hiding merging features, especially those like the ones 
we have in our sample that shine brightly in the mid-IR. 
We thus make them into our sixth, but least confident, 
merger category. 

The only objects that never enter our merger count are 
the spiral galaxies, since they are clearly isolated (that 
is non- merging) systems. 

We thus find that our sample spans the whole spec- 
trum of certainty, from secure mergers to faint ellipticals 
without any detected features. At the same time, we dis- 
cover that most of our sample lies in the middle of that 
range with signatures typical of, but not always exclu- 
sively associated with, a merger event. The most secure 
categories contain few of our merger candidates, and are 
thus highly incomplete. The most uncertain ones, on the 
other hand, show very little features and are thus fairly 
unreliable. In order to achieve the clearest possible pic- 
ture of the importance of merging in our sample, we need 
to attempt to strike a balance between completeness and 
reliabili ty. We discuss where that balance might lie in 
section |4.1[ together with more secure lower and upper 
limits. 
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Fig. 4. — Clear mergers (category 1). All of them are, by definition, in the early phases of merging: four are in phase III , and one is a 
triplet (see text for the definition of our morphological classes). The point source in the image of MIPS289 is a foreground star. 
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Fig. 5. — Highly probable mergers (category 2) separated by morphological class. For phases III & IV, both the image, and the residuals 
after subtraction of the main component are shown. Residuals for MIPS322 have been smoothed to better show its tidal features. This 
figure also serves to illustrate the range of morphological profiles displayed by objects of those two phases. 
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Fig. 6. — Likely mergers (category 3). Phase III objects in this category can have a higher luminosity ratio between their two components 
than their highly probable counterparts. They can also be separated by a larger distance, and less obviously connected. Phase IV objects 
possess tidal features characteristic of mergers, but they tend to be fainter than for those of the previous category, and can be sometimes 
incomplete. 
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but the shape of their features is not necessarily unique to the merging process. Phase V objects have very weak residuals, but these appear 
to form tidal features typical of mergers. Faint and compact objects have some of both. 



2.5.3. A word on automated classification techniques 

In recent years, automated galaxy classification using 
measur able parameters such as concentration and asym- 
metry ([Abraham et al. 1996; Cons eTice et al.|2003 ), Gini 
and M2U (|Lotz et al. 2004, 2008a L or other similaTcombi- 
nations (e .g. Scarlata et al. 20(T^ Law et al.|2007 Zamo- 
jski| [2008 ) have become the norm in high redshift galaxy 
surveys. Using these techniques, most mergers, in op- 
tically selected sam ples, can be identified by their high 
level of asymmetry ( Conselice et al.|2003[ ). In the case of 
ULIR Gs, the (?zm-M2(J scheme provides a cleaner sepa- 
ration ( Lotz et al.|2004 ), although the combination of the 
Gini coefficient and the asymmetry parameter appears to 
be, overall, the most sensitive to m erger signatures ( |Lotz| 



et al.|2004, 2008b Zamojski|2008|). 



Despite being undenia bly useful , especially for large 
samples, all of these automated classification techniques 
invariably miss a non-negligeable fraction of merger- 
induced ULIRGs, and fall short of the 99% mark inferred 



multiple nuclei systems are properly identified, but many 
of the singly-nucleated ULIRGs, i.e. the ones in more ad- 
vanced stages of merging, as well as mor e distantly sepa- 
rated pairs, fall outside the merger cuts (Lotz et al. 2004 



visually by Veilleux et al. (2002). In general, double or 



2008bj). This is due to the fact that in more advancec 
stages, merger signatures become increasingly faint com- 
pared to the remnant galaxy, while in distant pairs, the 
two objects tend to be treated separately, but not per- 
turbed enough yet, individually, to make the cut. In the 
case of advanced or late mergers, visual inspection, be- 
fore and after subtraction of a smooth profile, is usually 
the privileged, and perha ps the only effective, approach 
(e.g. | Veilleux et al.| ]20Q6). This is analogous to quasar 
hosts' lor w hich Fob -subtraction is crucial (e.g. Gabor 
et al.|2009 ). 

In section 3, we describe simulated observations we 
performed on local ULIRGs redshifted to z = 1. Ap- 
plying automated classification techniques to these red- 
shifted ULIRGs indicates that they become even harder 
to separate from the rest of the galaxy population by 
morphological parameters, with the best criterion (Gini- 
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— Possible mergers (category 5). They show some residuals, but one cannot tell for sure whether they originate from a merger 
not. Residual images have been slightly smoothed to help bring out the features. 



asymmetry) recovering only half of our 22 redshifted lo- 
cal ULIRGs listed in Table [I] In the rest-frame /-band, 
these sources have relatively smooth and broad tidal tails 
and bridges. When put at high redshift, because of sur- 
face brightness dimming, these tend to get buried under 
the noise . Morphological parameters then become in- 
creasingly derived from the central components which, 
themselves, appear smoother due to the decreased reso- 
lution. These effects tend to draw those objects towards 
the same morphological space as that spanned by nor- 
mal disks and bulges, therefore making mergers harder 
to identify. This is what happens to three of the four 
sources illustrated in Figure [9] (IRAS 22491-1808 being 
the only one to retain merger-type parameters at high 
redshift), and this despite the fact that they still show 
good indications of being mergers, even at z = 1. 

The same processes appear to affect our observed sam- 
ple. In addition, a large fraction of our objects possess 
prominent bulges whose signal overwhelms that of merg- 
ing features, even in bright objects. This is then re- 
flected in the value of their morphological parameters, 
which overlap considerably with that of normal galaxies 
and often fall short of the merger criteria. These objects 
are often at lower redshift (z < 1) where they are seen 
in their rest-frame near infrared. As a result, we find 
many clear mergers not making any of the morphologi- 
cal cuts discussed above, and, typically, only around one 
third of our sample would appear as mergers using these 
techniques. We find that those objects that do make the 
cuts are really only the most obvious one: mainly those 
with double or multiple nuclei. This is extremely limita- 
tive. We would rather include, in our study, objects in 
all stages of the merging process. We therefore choose to 
forego, in this paper, the use of those parameters on our 
24/im-selected sample as we feel they do not reflect well 



its structural richness. 

The situation is very different for our control sample. 
The latter is composed primarily of disky galaxies, and 
resembles other optically /NIR-selected sample for which 
autom ated techniques such the CAS system (Conselice 
|2QQ3[ ) were designed. It is, therefore, to no surprise that 
we find a much better agreement between visual and 
automated classifications with that sample. A simple 
asymmetry cut suggests a merger fraction of 20% in our 
control sample while a combination of Gini and asym- 
metry yields a merger fraction of 28%. For comparison, 
the merger fraction inferred visually is 25 % when count- 
ing up to confidence level 4 (section 4.1 Fig. 10). Both 



estimates are domin ated by objects with H ?5l8-20 as 
explained in section |4.1| This suggests that our visual 
classification is indeed reliable up to confidence level 4, 
so that we can feel confident using it in the regime where 
automated techniques appear to break down, that is the 
one to which most of our 24/im-selected sources belong. 

Finally, we note that the merger fraction in our con- 
trol sample is higher than most studies would suggest 
for a population at z < 1 flLotz et al. |2008a Kartal- 
tepe et aLl|2007| Cassata ~et al.||2005| |Conselice et al. 
2003). BasecTon the local J-band luminosity function 



([(Jole et al. 2001), z = 0.15-1 is where we expect to 
find mos t of our H = 18-20 sources. We showed in sec- 
tion [2^21 however, that we do have an excess of sources in 



our images (that is without counting the 24/im-sources) 
at these magnitudes compared to the average i/-band 
number counts. This excess must come from the envi- 
ronment these 24/im-sources live in, and we therefore 
attribute our higher merger fraction to that over-density 
of galaxies around 24/im sources. 

2.6. Mid-IR Spectral Diagnostics 
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The mid-infrared part of the integrated spectrum of 
a galaxy contains imprints of the physical processes 
at the origin of its infrared luminosity. These in- 
clude PAH emission complexes and silicate absorp- 
tion/emission bands, as well as continuum emission from 
hot (T ~ 1000/T) dust, if present. We discuss below how 
these features inform us on the nature of our sources and 
how we measure them. 

2.6.1. PAH Equivalent Widths 

Polycyclic aromatic hydrocarbon (PAH) molecules ex- 
ist and are excited primarily in photo-dissociation regions 
around OB associations. This makes them good tracers 
of star formation. When excited, they reradiate their en- 
ergy in a number of narrow bands at wavelengths mainly 
between 6.2 and 12.7 jam. These, in turn, straddle the 
broad peak of hot dust continuum emission that spans 
the entire mid-IR range. Since, in external galaxies, the 
primary source of dust heating is, when present, the ac- 
tive galactic nucleus, the ratio of the PAH flux to that 
of the underlying continuum, or equivalent width, gives 
us a proxy for the relative importance of star formation 
versus AGN activity. 

In this paper, we use both the equivalent width of the 
7. 7 jam and 11.3/xm PAH features as indicators of the re- 
spective contribution of the starburst and AGN compo- 
nents to the mid-infrared flux of our objects. We perform 
our measurements using the M arkov Chain Monte C arlo 
fitting procedure d escribed in Sajina et al. (2007) and 
|Sajina et aLlfl2006| ). 

Comparison of equivalent width measurements from 
our method with full SED-fitting indicates that objects 
with an EW (11. 3 fim) > 0.8 or an EW(7.7jim) > 1.2 are 
powered at > 80% by star formation, whereas those with 
an EW (11. 3 fim) < 0.1 or an EW(7.7jim) < 0.15 have 
a < 20% contribution from star formation. Through- 
out this paper, we will refer to the former category 
of sources as starburst-dominated, and to the latter as 
AGN-dominated. Intermediate objects are given the ep- 
ithet "composites". 

Lastly, we mention that our boundary between AGN 
and composite systems is slightly blurred due to the fact 
that our upper limits of detectability of PAH features 
often lie around or just above the division line between 
these two categories. In this paper, we choose to con- 
sider all non-detections as AGN-dominated systems. We 
recognize that a few of them might actually be low-end 
composites, but none of our conclusions rely on knowl- 
edge of the exact ratio of the number of AGNs to low-end 
composites. 

2.6.2. The 9.7 jam silicate feature 

Silicate grains are ubiquitous throughout the interstel- 
lar medium. They show a broad emission profile centered 
at 9.7 jam when heated to temperatures of a few hundred 
kelvins or more, but the profile is often rather seen in ab- 
sorption as colder grains screen the line-of-sight towards 
a hotter source. In distant galaxies, this hot source is the 
galaxy's active n ucleus. In the scenar io of quasar forma- 
tion proposed by Sanders et al. (19881), the latter are born 
in the dense, highly obscured cores of coalescing galaxies. 
As they l ight up, they start pushing away surrou nding 
material (Silk & Reesl 19981 |Di Matteo et al.||2005D, and 



(Jonsson et al. 2006). The presence of deep silicate ab- 
sorption can therefore be indicative of this early dense 
phase. 

In this paper, we use the optical depth of the 9.7 fim 
silicate feature, rg.7^ m , primarily to look for these young 
obscured AGNs in our sample, but also as a general proxy 
for the level of obscuration in our objects. It is defined 
as: 



7"9.7/xra 



In 



fcontX9'7fim) 
/(9.7/xm) 



(1) 



We use the method of iSajina et all ([20071 to fit the 
optical depth of the silicate absorption feature simul- 
taneously with the PAH and continuum fluxes used to 
obtain PAH equivalent widths. This procedure assumes 
a screen geometry. The optical depth at 9.7 jam relates 
to t he silicate strength often in use in the literature 
(e.g. Spoon et al. 2007| ) through the following relation: 

7"9.7/xm = —1.4 X S s u. 

3. SIMULATIONS 

A key task of any high redshift study is to establish a 
comparison with more nearby samples. Because merg- 
ers and merger remnants are such defining trademarks 
of local ULIRGs, it is imperative, and the primary focus 
of this paper, to ask whether these traits are as ubiq- 
uitous among their high redshift counterparts. Surface 
brightness dimming and lower resolution, however, af- 
fect our ability to detect merger signatures and identify 
mergers at high redshift. We address these issues in this 
section using simulated HST/NICMOS observations of 
local ULIRGs redshifted to redshifts of one and above. 

We artificially redshifted a total of 22 local ULIRGs 
and merge rs that we selected primarily from the IRA S 1- 



Jy sample jKim fc Sanders 1998; Veille ux et aL]2002 | ) and 
the sample of Surace et al. (1998, 2000) in such a way as 
to cover a wide range of merger configurations and mor- 
phological types, as well as both warm (/25 M m//60 M m > 
0.2) and cold (f25fim/ feOfim < 0.2) mid-IR colors. We 
used mainly archived HST/WFPC2 F814W data when 
available or, other wise, /-band ground-based images 



levels of obscuration are, consequently expected to drop 



from Surace (1998]) taken with the UH2.2m telescope. 
Those filters were chosen to match the rest-frame band 
observed by the NICMOS/F160W filter at z = 1. For 
HST images, we used the Tiny Tim software to model 
the PSF, whereas for ground-based images, we used stars 
directly available from the image itself. The list of galax- 
ies we simulated can be found in table T] along with the 
provenance of the data used for each object, basic infor- 
mation about the object such as its morphological class, 
redshift, magnitude and mid-IR color, and its classifica- 
tion after redshifting. Figure [9] shows examples of simu- 
lated images of local phase III and phase IV ULIRGs. A 
full description of our redshifting procedure can be found 
in Appendix A. 

We elected to simulate most of our objects at z — 1 
for two reasons: 1) it corresponds to thepeak of the red- 
shift distribution of our sample (Figure |2|, and 2) /-band 
images, which correspond to the rest-frame wavelengths 
we observe at that redshift in our data, were the ones 
most easily available to us. Although we also have a sig- 
nificant number of objects, in our sample, at z ~ 2, we 
actually find that the merger nature of those objects is 
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Fig. 9. — Examples of simulated observations of redshifted local ULIRGs. The color scale is inverted in residual images, and they have 
been smoothed to help bring out faint features. 
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TABLE 1 
Redshifted objects 



Object Name a Instrument Filter Magnitude Redshift b Absolute mid-IR Simulated Classification after 

Magnitude color c Redshift(s) redshifting d 



Phase III objects 



IRAS 08572+3915 


HST/WFPC2 


F814W 


15.59 


0.058 


-21.55 


warm 


1 


III (2) 


IRAS 14348-1447 


UH2.2m 


I 


14.58 


0.082 


-23.36 


cold 


1 


III (1) 


IRAS 03521+0028 


HST/WFPC2 


F814W 


17.49 


0.152 


-21.96 


cold 


1 


Irregular 6 


Mrk 463 


HST/WFPC2 


F814W 


13.34 


0.0504 


-23.46 


warm 


1 


III (2) 


IRAS 12112+0305 


UH2.2m 


I 


15.05 


0.072 


-22.59 


cold 


1 


III (1) 


IRAS 22491-1808 


UH2.2m 


I 


15.00 


0.078 


-22.82 


cold 


1, 1.7 


III (1) 


IRAS 23498+2423 


HST/WFPC2 


F814W 


16.81 


0.212 


-23.50 


cold 


1 


III (3) 



Phase IVa objects 



Arp 220 


HST/WFPC2 


F814W 


12.25 f 0.0185 


-22.28 


cold 


1 


F&C§ (5) 


IRAS 15206+3342 


HST/WFPC2 


F814W 


15.61 0.125 


-23.22 


warm 


1 


IV (3) 


IRAS 02021-2104 


HST/WFPC2 


F814W 


16.29 0.115 


-22.48 


warm 


1 


IV (4) 


IRAS 09039+0503 


HST/WFPC2 


F814W 


16.27 0.126 


-22.59 


cold 


1 


IV (3) 



Phase IVb objects 



Mrk 273 (UGC8696) HST/WFPC2 F814W 13.16 f 0.0382 -22.97 cold 1 IV (2) 

IRAS 15250+3609 UH2.2m I 15.19 0.0555 -21.78 cold 1 V (5) 

UGC 5101 UH2.2m I 14.60 0.0394 -21.60 cold 1 F&CS (6) 

Mrk 231 HST/WFPC2 F814W 12.49 0.0426 -23.92 warm 1 IV (4) 

Mrk 1014 HST/WFPC2 F814W 14.84 0.162 -24.44 warm 1 IV (3) 

IRAS 23365+3604 HST/WFPC2 F814W 14.48 0.0634 -22.79 cold 1 F&CS (5) 

IRAS 05189-2524 HST/WFPC2 F814W 13.91 0.0427 -22.51 warm 1 IV (4) 



Phase V objects 



IRAS 01003-2238 
IRAS 07598+6508 
IRAS 13218+0552 

Mrk 771 



HST/WFPC2 F814W 17.15 0.117 -21.61 warm 1 Regular bulge 

HST/WFPC2 F702W 14.40 0.149 -24.84 warm 1.32 Regular bulge 

HST/WFPC2 F702W 17.24 0.206 -22.55 warm 1.32, 1.8 Regular bulge 

HST/WFPC2 F606W 14.95 0.064 -22.36 warm 1.67 V (5) 



a Objects in bold are part of the IRAS 1-Jy sample 
b Corrected to the rest-frame of the CMB 

c Objects with a ratio of /25/^m//60/^m > 0.2 are considered warm. Otherwise, they are said to be cold. 
d Confidence levels (see §2.5.2) are shown in parentheses 

e IRAS 03521+0028, when redshifted, looks nothing like any of our objects. It rather looks like a diffuse cloud, a scaled up version of 
the LMC. We therefore call it irregular, but would not have associate such morphology with a merger had we encountered it in our 
sample. 
f From ! 
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Surace 



6z, (Jompact 



(1998) with the UH2.2m telescope. 
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often clearer than that of the z ~ 1 sources. This is most 
likely due to the fact that, at z > 1.5, objects in our 
sample have infr ared luminosities in excess of 1O 12 ' 5 L 



(see section 4.4), and are, therefore, more likely to be 
more major mergers. Furthermore, we show later that, 
past z = 1.5, local ULIRGs are no longer good analogs 
of our bright 24/im-sources, whether morphologically or 
in terms of their spectral properties. It is, therefore, un- 
clear whether local ULIRGs can be used as calibrators 
for our highest redshift sources even if we were to per- 
form simulated observations of them at z = 2. For these 
reasons, we feel that z = 1 is the most relevant redshift 
to put these objects at, and that this provides us with 
adequate and sufficient information for our purposes. 
The results of these simulations reveal that most phase 

III mergers are still identifiable as such at high redshift s. 
They demonstrate the same to be true for two thirds of 
our simulated phase IV objects as well. Table [I] also pro- 
vides evidence, however, that our phase IV is narrower 
than its local homologue, and that, at least some, phase 

IV objects are likely to, instead, fall in the phase V cat- 
egory in our sample. 

It also emerges from our simulations that objects we 
classify as Faint & Compact at high redshift are always 
identified as phase IV mergers at higher resolution and 
higher signal-to-noise. These objects typically possess a 
dense core, and are surrounded by faint and diffuse tidal 
structures. At high redshifts, only the core is seen, which 
is what gives them their "faint & compact" appearance. 
It is this observation that motivates our decision to asso- 
ciate these faint & compact objects with the coalescence 
phase, together with the advanced (phase IV) mergers. 
By doing so, we are further able to recover an even 
greater fraction of these mergers. 

The larger difficulty arises with older, phase V, merg- 
ers. We simulated high-redshift observations of four of 
such objects and, indeed, find that, in most cases (three 
out of four), the merger signature is lost in the process. 
This demonstrates that our data is not very sensitive 
to the latest stages of the merging sequence and that 
late mergers could masquerade as regular elliptical galax- 
ies. When signatures are detected, they further consist 
mostly of lopsidedness or asymmetries which are diffi- 
cult to link uniquely to a merger event. We take all this 
into account in our estimate of the merger fraction in 
section 14. 11 

Although our simulations have shown that, for the 
most part, local ULIRGs can still be identified as merg- 
ers at redshifts of z > 1, they do not address the more 
difficult question as to whether, or to what extent, non- 
mergers could falsely pose as mergers at lower resolution 
and signal-to-noise. It is, in fact, possible to imagine an 
alternative scenario for nearly all of our merger-identified 
objects, whether it is chance superposition accompanied 
by coincidentally misleading distortions giving the im- 
pression of an early merger, or fly-by interactions giv- 
ing rise to tidal tails and posing as an advanced merger, 
or yet the presence of an undetected dwarf companion 
causing the primary object to display a lopsided pro- 
file and asymmetry characteristic of phase V objects, are 
all within the realm of possibility. We argue, however, 
that such serendipitous configurations, at least among 
phases III and IV, must occur far less frequently than ac- 
tual mergers, most particularly in bright 24/ira-selected 



sources, and that a full quantification of the level of con- 
tamination is, therefore, unnecessary at this stage. In- 
stead, we chose to assign a confidence level to each of our 
objects and quote our merger fraction as a function of 
confidence level rather than in absolute terms. This will 
allow interpretation of our results to easily carry into the 
future as calibration of our confidence scheme, in terms 
of completeness and reliability, clarifies itself, and the 
uncertainty on the merger fraction narrows down. We 



describe this approach in more details in section 4.1 



As a final note, we turn our attention once again to our 
Faint & Compact objects, and mention that even though 
isolated disks that are faint and compact are known to 
exist at high redshifts ( Guzman et al.|1 998; Noeske et al. 
|20Q6[ ), they tend to be very blue ( |Phillrps et al.|l997| ), im- 
plying that their faintness is intrinsic and not due to ex- 
treme extinction. They would, therefore, not make it into 
our bright 24/im flux-limited sample. The mid-infrared 
properties of our faint & compact objects presented in 
section 4.5 and 4.6 would, further, be very difficult to 
reconcile with isolated star-forming disks. This suggests 
that compact isolated disks are distinct from our Faint & 
Compact sources and that contamination is not likely to 
be a major issue. We nevertheless, treat these objects as 
mergers with a lower level of confidence because, despite 
the fact that all the evidence seems to paint a consistent 
picture, it still remains circumstantial in nature. 

4. RESULTS 

4.1. Merger fraction 

To describe the merger fraction of a sample in abso- 
lute terms is a perilous proposition. It is, however, the 
first question that arises from a data set such as ours 
and, therefore, needs to be addressed. We know that, 
locally, interacting galaxies are readily identified visu- 
ally and that high resolution observations can still reveal 
telltale signatur es of past events even at late stages o f 
the process (e.g. |Surace et al.|1998 Veilleux et al.|2 QQ6). 
This confidence, however, decreases exponentially with 
redshift as merging signatures rapidly blur and fade, and 
calibration of lower resolution and signal-to-noise obser- 
vations against the multitude of galaxy morphologies and 
configurations observed in the nea rby and distant Uni- 
verse, despite nume rous efforts (e.g. Conselice 2003; Lotz 
|et al.||2004} |2008b[ see also sections~ [2!5.3| and 3), re- 
mains, toaate, incomplete. Consequently, one has to 
rely on visual classification where opinions might differ. 
In an effort to be transparent, we decided to approach the 
questio n in st eps by splitting our sample, as discussed in 
section 2.5.2, into confidence classes, and then account- 
ing for the merger fraction cumulatively from the most 
to the least certain of them. The results are shown in 
Figure [To] We provide our own interpretation of these 
results below, but our method also allows them to be 
more easily re-interpreted in light of improved calibra- 
tion constraints in the future. 



Figure [10] illustrates our full accounting of the merger 
fraction by confidence level and morphological class. It 
also shows a similar curve for the 1- Jy sample as if it were 
observed at high redshift with the same kind of data, as 
well as one for our control sample. The curve for the 
redshifted 1- Jy sample was produced using the results of 
our simulations listed in Table [I] applied to the intrinsic 
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Fig. 10. — Cumulative merger fraction from the most to the least 
certain confidence class. Classes 1 through 5 are illustrated with 
examples in Figures [4] to [8] Class 6 are regular, but faint galaxies, 
and could therefore oe hiding features. Our best estimate of the 
merger fraction is 80% (see text for details). The grey area repre- 
sents the range of merger fractions obtained for our control sample 
from automated classification techniques, illustrating that our vi- 
sual classification is reliable up to confidence level 4. Automated 
c lassific ation techniques fail for ULIRGs at high-z as discussed in 
§ |2.5.3| and hence comparisons are not shown for either ou sample 
or the redshifted 1-Jy sample. 



morphological composition of the sample. Because the 
magnitude distribution of objects in our control sample 
is wildly different from that of our 24/im-sources (cf. Fig- 
ure]^]), a direct comparison of the merger fraction in those 
two samples would not be meaningful. Instead, to derive 
the merger fraction shown in Figure [10} we weighted the 
contribution of every object in our control sample by the 
ratio of the number of objects within that same magni- 
tude bin in our 24/im-selected sample to its number in 
our control sample. As a result, most of the weight, in 
the merger fraction of the random sample lies in objects 
with magnitudes ranging from H ~ 18-20. We use these 
two curves to inform our estimate of the merger fraction 
in our 24/im-selected sample. For example, a compari- 
son of the cumulative merger fraction curve to that of 
the redshifted 1-Jy sample immediately reveals a strik- 
ing similarity, from which we can already conclude that 
our sample is at least consistent with a high merger frac- 
tion. More importantly, however, a high merger fraction 
appears to be necessary to explain the gap between our 
24/im-selected and control samples, the latter of which 
we find to be al ready composed of ~ 25% mergers by all 
measures (see § 2.5.3). We argue below that, in fact, the 
merger fraction among our bright 24/im-selected galaxies 
is unlikely to lie outside the 62-91% range. 

As a start, we can set a reasonable lower limit on our 
merger fraction simply by using the fact that, as demon- 
strated in Figure 10 our 24/im-selected sample bears 
much closer resemblance to the redshifted 1-Jy sample 
than to the control sample. We can, therefore, be confi- 
dent that the merger fraction in our sample lies beyond 
the mid-point between those other two samples, that is 
above 62%. 

In order to arrive at a more accurate estimate of the 



true merger fraction, however, we need to rely a little 
bit more on our intuition to find a balance between com- 
pleteness and reliability. This question remains, there- 
fore, somewhat subjective and one can, naturally, form 
one's own opinion based on the images shown in Fig- 
ures]^ to [8] However, one still needs to explain the ob- 
served morphologies of our objects, and trying to do so 
without invoking a merger, we more often than not found 
ourselves facing much more improbable scenarios. 

We have two pieces of evidence that allow us to b e a lit- 
tle more quantitative. First, as discussed in section[2.5.3 



there is excellent agreement between our visual classi- 
fication and automated classification techniques when 
counting up to confidence level 4, in our control sample. 
(Unfortunately, because of the limitations of those tech- 
niques, this test is not applicable to our 24/im-selected 
sample). Then, an examination of the ACS and WFPC2 
images of those 24/im-selected objects that fall within 
the area covered by those instruments, confirmed, that 
is removed any doubt as to, the merging nature of many 
of these objects. It even revealed the merging nature 
of two sources that we had thought, previously, to be 
disturbed disk galaxies (MIPS168 and MIPS8328). On 
the other hand, it never reversed any of our NICMOS 
merger identifications. This unidirectionality, along with 
the agreement with automated techniques in our con- 
trol sample and the difficulty of explaining our observed 
morphologies through alternate means, all suggest that 
confidence levels one to four, although increasing in un- 
certainty, must still be primarily populated by mergers, 
and that contamination probably amounts to no more 
than a handful of objects. We also argue, based on our 
redshifting analysis of the 1-Jy sample, that many, and 
possibly the majority, of category 5 and 6, 24/im-selected 
objects also have their origin in an ongoing or recent 
merger event, and that the number of mergers in those 
two categories, thus, likely exceeds the number of false 
positives in the lower levels. We therefore arrive at the 
conclusion that the correct balance between complete- 
ness and reliability should lie somewhere in the vicinity 
of confide nce le vel 5, and that is where we draw the line 
in Figure lCf^ 



Our best estimate of the merger fraction thus lies at 
%. We can also arrive at the same number by making 
the simplifying assumption that all confidence class 1 
through 4 objects are mergers, while only half of the 
class 5 and 6 ones are. Since all of the phase I through 
IV mergers belong to confidence classes 1 through 4 and 
almost all of the faint & compact and late mergers belong 
to classes 5 and 6, accounting for all of the former group 
and half of the latter also yields a similar fraction. 

In order to get a handle on the uncertainty associated 
with this number, it is useful to examine what our merger 
fraction would be under different assumptions. For ex- 
ample, if we were to draw the line at category 4 and ex- 
clude all category 5 and 6 objects from the merger count, 

11 We note that the argument about confidence levels 5 and 6 
objects does not hold for our control sample as most objects in that 
sample are extended disk galaxies for which we have no compelling 
reason to think that the presence of residuals should come from 
a merger more than from any normal asymmetry or structure in 
the disk. The correct limit for that sample should, therefore, be 
drawn at level 4 and no further. The merger fraction at that point 
reaches 25%. 
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we would obtain a merger fraction of 72%. Counting all 
category 1 to 3 objects, but only half of those in cate- 
gories 4 to 6, would similarly yield a 69% fraction. Such 
assumptions, although seemingly less likely, are still re- 
alistic, and cannot, therefore, be ruled out. If we were 
to draw the line even further down, somewhere between 
categories 3 and 4, which amounts to counting all ob- 
jects in confidence levels 1-3 plus about half of those in 
level 4, our inferred merger fraction would be ~ 60%. 
60% is also the fraction one would get by assuming a 
20% contamination in the first four categories while still 
rejecting all category 5 and 6 objects. Such skepticism 
would imply that a large number of objects are simply 
masquerading as mergers, and that 40% of our observed 
morphologies would have to be explained through alter- 
nate means. This becomes very hard to justify, especially 
in the view of the fact that ~ 60% of our objects possess 
an important bulge component. We feel that such a point 
of view is very much at the limit of what is reasonable to 
assume. We thus find that our earlier lower limit of 62% 
on the merger fraction derived from a simple comparison 
of our sample with the redshifted 1-Jy and control sam- 
ples corresponds well to what we intuitively infer from 
direct inspection. We therefore retain that number as 
our lower limit. 



At the other end, Figure 10 demonstrates that a sample 
of nearly 100% merger-driven ULIRGs, namely the 1-Jy 
sample, can have a cumulative merger fraction curve very 
close to that of our sample. It is thus not impossible that 
all of our merger candidates actually be mergers. For this 
reason, we use the maximum value of the curve for all 



mergers in Figure 10 as our upper limit. It reads 91%, 
and excludes only isolated spiral galaxies. The range 
of reasonable values for the fraction of mergers in our 
sample is, thus, 62 to 91%. 

Historically, infrared-luminous galaxies have been di- 
vided between LIRGs and ULIRGs. For comparison pur- 
poses, we thus also performed separate accountings for 
those two categories individually using the same dividing 
lines. The numbers add up to a 57-80% merger fraction 
among LIRGs with a best estimate of 71%, and to a 
merger fraction between 65% and 96% among ULIRGs 
with a best estimate at 87%. The reason why our merger 
fraction is so large among LIRGs is that most of our 
LIRGs have luminosities just below 1O 12 L , the mean 
luminosity being (log Lir/Lq) = 11.63. The numbers 
we find for ULIRGs, on the other hand, indicate that, 
like their low redshift counterparts, they too have their 
origin primarily in mergers. 

4.2. Interlude: postulate for results to follow 

In the sections that follow, we look for correlations be- 
tween various properties of our galaxies and their mor- 
phological class. Incorporating confidence levels into 
all of our relations would be impractical. We there- 
fore choose to leave them out and include in our sub- 
sequent analysis all objects of every morphological class 
irrespective of their confidence level. We argue this is 
justified, because the trends discussed in this paper are 
robust and would not be affected by the reshuffling of 
a few objects around morphological classes. They fo- 
cus mainly on the distinction between early and more 
advanced mergers which are easily distinguishable by 
whether the source consists of two objects/nuclei or one. 



Cross-contamination between these two classes of objects 
is therefore likely to be minimal. Moreover, Figure 10 
demonstrates that early mergers and advanced mergers 
have a very similar curves as for their cumulative merger 
fraction as a function of confidence level, implying that 
eventual interlopers are likely to affect both categories 
equally (though, as argued in the last section, their num- 
ber is expected to be very small). What's more, the 
objects associated with the least confidence to a merger 
event (faint & compact objects, and late mergers) have 
very uniform spectral properties, so that the conclusions 
remain the same for those categories irrespective of which 
objects one chooses to include or not. 

This assumption, or approximation, we choose to pur- 
sue greatly simplifies later analysis, and allows us to re- 
tain the focus on the more physical merging sequence. 

4.3. Morphological split along the merging sequence 

Having established the overall merger fraction, we now 
look into the distribution of our mergers along the merg- 
ing sequence (see section 2.5.1 for a description of our 
classification). This distribution is especially meaningful 
when laid in comparison with other samples. For that 
purpose, we continue to u se the 1-Jy morpholog ical sam- 
ple of local ULIRGs from Veilleux et al. (2002), and the 
results of our simulations to compare with local galaxies 
as if they were observed at redshift one with the same 
kind of data. 

We show in Table [2] the morphological split of our sam- 
ple (also represented in graphical form in Figure 10 ) and 
compare it to that of the 1-Jy sample. We list num- 
bers for both the intrinsic and redshifted 1-Jy sample, 
the percentages for the latter being obtained by applying 
the proportions found in table [I] to the sample's intrinsic 
morphological distribution. Finally, we divide our sam- 
ple into two redshift bins separated at z = 1.5 and show 
the morphological composition of each sub-sample. Our 
full sample roughly splits into 1/3 early mergers, 1/3 ad- 
vanced mergers and 1/3 other classes. Among the last 
third, the most represented types are old mergers (14%), 
isolated spirals (9%) and faint & compact objects (8%). 

Comparison between our objects and local ULIRGs 
reveals that not only both samples are composed, to 
a high fraction, of mergers, but that all phases of 
the merging process are represented in relatively simi- 
lar proportions. In particular, both are dominated by 
singly-nucleated merger remnants (coalescence and post- 
coalescence mergers). We find, however, a number of 
subtle differences between the two samples, and address 
these differences below. 

The main difference is in the number of isolated spirals. 
Whereas all but one of the 1-Jy galaxies are mergers, we 
have, in our sample, 12 isolated disks (7 face-on and 5 
edge-on). We also have 11 faint & compact objects (8 of 
which show residuals after subtraction of the main profile 
and 3 that do not). Some of these objects could also be 
intrinsically isolated (i.e. not merger remnants). This 
difference in the number of isolated disks is explicable 
by the fact that some of our lowest redshift objects are 
sub-ULIRG (see Figure [IT ) . 



The largest discrepancy among merger classes occurs 
in phase V objects, for which we would have expected, 
from our redshift ing of local ULIRGs, to find a lower frac- 
tion, and to observe, instead, many more objects with no 
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TABLE 2 

Morphological content of different samples & 



Merging sequence 


Morphological class 


1-Jy sample 


Redshifted 


This 


z < 1.5 


z > 1.5 
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1-Jy sample 


sample 


only 


only 
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1% 
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Isolated objects 


Edge-on disks 
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6% (4) 





1 




Irregular galaxies 




4% 
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3% (4) 


4% (3) 


2% (1) 


o 




First contact (phase II) 








2% (3) 


3% (2) 


2% (1) 





Early mergers 


Pre-mergers (phase III) 


39% (46) 


35% 


27% (36) 


12% (8) 


46% (25) 


3 




Triplets 


4% (5) 


4% 


3% (4) 


3% (2) 


4% (2) 







Advanced 














Coalescence 


Mergers (phase IV) 


44% (52) 


29% 


31% (41) 


37% (25) 


28% (15) 


1 




Faint & Compact 




12% 


8% (11) 


9% (6) 


7% (4) 


1 




Old mergers (phase V) 


12% (14) 


7% 


14% (18) 


12% (8) 


9% (5) 


5 


Late mergers 


Regular bulges 





9% 


3% (4) 


4% (3) 


2% (1) 







Pure point sources 







1% (1) 








1 



a Numbers without percentages represent total number of objects. Ellipses indicate that morphological class is not used for 
that sample. 



signs at all of interaction (that is regular bulges). How- 
ever, we see quite the opposite, with phase V objects be- 
ing far more common than simple bulges. This suggests 
that our objects carry stronger merger signatures than 
do local ULIRGs. We see that effect among advanced 
mergers as well when visually comparingtidal tails and 
streams of objects in our sample (Figures|5]toj7|) to those 
of the redshifted 1-Jy objects (Figure [9]). We find that 
those of our sample tend to be less diffuse, higher surface 
brightness, and readily detectable. Among our redshifted 
phase IV local ULIRGs, on the other hand, only Mrk273 
still shows a clear tidal tail at z = 1. More typical galax- 
ies, such as those shown in Figure [9j display only very 
faint streams and tails that are hard to detect but for 
smoothing of the residuals. We speculate that this effect 
is one more manifestation of the enhanced gas fractions 
and star formation rates that exist at those redshifts. 

Although, overall, our sample shows a number of early 
mergers very close to that of the redshifted 1-Jy sample, 
the proportion of these objects is three times as large 
at high redshift than it is at z < 1.5. We defer discus- 
sion of the large number of early mergers at high redshift 
until the next section, where we argue for a redshift evo- 
lution in the morphology of ULIRGs. The low number 
of early mergers among our z ~ 1 objects compared to 
that of the 1-Jy sample, on the other hand, can proba- 
bly be explained by the fact that the two samples were 
selected at different rest-frame wavelengths, namely at 
A ~ Yl\im for our objects and at 60/im for the 1-Jy 
sample. Many objects in the 1-Jy sample have fairly 
low mid-IR fluxes (those with high / (60/im) //(25/ira) 
ratios). Such objects are less likely to make it into our 
sample. They are, however, predominantly early-stage 
mergers flVeilleux et aL|2002| ). 



4.4. Redshift evolution of morphology 



In Figure 11 we show the total infrared (3-1000/im) lu- 
minosity as aiunction of redshift for all galaxies of known 
redshift and luminosity in our sample. The distribution 
is typical of a flux-limited sample in which the more lumi- 
nous sources are at higher redshift and vice versa. The 



sample seems to split naturally, though, at z = 1.5 or 
Lir = 1O 12 5 L , with clearly different populations on 
each side. We explore the differences between those two 
populations, as well as with local ULIRGs, further in this 
section. 
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Fig. 11. — Distribution of IR-luminosity for all objects of known 
redshift in our sample as a function of redshift. The average 3cr 
uncertainty is shown in the upper left corner. The light gray line 
represents the I R-luminosity at which star formation is predicted, 
from models of |Hopkins et ah] ( |2010| ), to transition from quies- 
cent to merg er-dr ivem (Jomparison between data and models is 
discussed in 8 15.51 



Figure [12] represents the morphological distribution of 
objects in our sample, split in high and low-redshift bins 
at z = 1.5, with the addition of a no-z column for objects 
with unknown redshift. It shows a clear shift from singly- 
nucleated mergers at low redshifts to earlier phase merg- 
ers at z > 1.5. This is true even if all objects of unknown 
redshift were at high-z, as illustrated in the last column. 
Splitting the sample in luminosity at Lir = 10 12 ' 5 L© 
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instead of in redshift yields identical results. 
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Fig. 12. — Number of objects of each morphological class in two 
redshift bins. Also distribution of morphologies for objects with- 
out known redshift. There are more early stage mergers at high- 
redshift even if all objects without redshifts were put in the high-z 
bin. 



So are we seeing redshift evolution or a luminosity se- 
lection effect? We argue that these observations cannot 
be explained by luminosity, since they go opposite to the 
luminosity trend known locally, whereby more luminous 
objects tend to be in more a dvanced stages of merging 



(Veilleux et al. 20021 |2009ab. Neither can they be ex- 



plained by an increase in the importance of star forma- 
tion activity since we actually find it to decrease with 
redshift in our sample, as illustrated in Figure [13] which 
demonstrates the smaller PAH equivalent widths of our 
high redshift objects compared to that of our lower red- 
shift objects; this, despite the fact that, at z ~ 2, sources 
are selected at the rest-frame ~ 8/ira where the presence 
of PAH features could potentially significantly increase 
their observed fluJ^l From the relation between redshift 



and luminosity in our sample (Figure 11), we can, inci- 
dentally, deduce that our objects follow the same trend 
as found in studies of local ULIRGs where PAH EWs 
have been shown to decrease, on a verage, with luminosity 



11U V kyv^V^ll U11UYV 11 UV U.W1 VUJIJV^ , VJ11 UJVV/lLl^V^ V V 1U11 1 UllllllVUl KJ J 

([Lutz et al.|1998||Tran et al |2001l [Veilleux et al.|2Q09b| |. 

The analogy with local ULlKGs, however, would predict 



would predict 



local ULIRGs, 
that at higher luminosities and lower EWs, which, m our 
sample, occur at higher redshifts, we find more advanced 
mergers. Since we observe exactly the opposite, we must 
conclude that there is intrinsic morphological evolution 
with redshift (which is unlike PAH EWs whose relation 
with redshift seems to be driven primarily by luminos- 
ity). This morphological shift towards early mergers at 
z > 1.5 therefore implies that, at those redshifts, more 
activity, and especially more AGN activity as testified to 
by the low PAH EWs of those objects, occurs in early 
stages of merging. We explore possible causes for this 
shift in the discussion section. 

13 The presence of 11.3/ira and 12.7/ira PAH features can also 
favor selection among lower redshift objects, but those lines typ- 



4.5. Relation between mid-IR properties and 
morphology 

In this section we look at the relation between the 
physical origin of the mid-IR flux of our objects, that 
is whether it is from star formation or AGN activity (or 
a combination of both) , and their observed morphology. 
As mentioned in section 2, we use both the EW of the 
11.3/ira and 7.7/ira PAH features as indicators of the re- 
spective contribution of the two phenomena. 

In Figures 14 and[l5]we show the distribution of equiv- 
alent widths split by morphological class for our z < 1.5 
and z > 1.5 samples respectively. Figure [14] demon- 
strates an evolution in the median PAH EW along the 
merging sequence, from high PAH EWs in isolated spi- 
rals and early mergers to low EWs or no PAHs at all 
in old mergers and regular ellipticals. This indicates a 
progression from star formation to black hole accretion 
activity as the merging process advances, much like the 
one o bserved in local ULIRGs (Figure 16 Veilleux et al. 
2009a). Just as in local ULIRGs, however, this evolution 
is statistical in nature: all sources do not go through that 
exact sequence. On the contrary, we observe a broad 
distribution of equivalent widths in most morphological 
classes, indicating a large variability. We find the same 
trends in Figure [l5j albeit on a narrower span of the 
merging sequence. 

Comparison of Figures [l4j [15] and 16 reveals that 
the PAH equivalent widths oitne three samples are 
overall shifted with respect to one another. This 
shift, however, is mainly a reflection of the fact that 
these samples cover different luminosity ranges, from 
logL/^/L© ~ 11.5-12.5 in our low redshift sample (high- 
est EWs) to log L IR /L Q ~ 12-12.5 among 1-Jy ULIRGs 
to \ogLm/ L@ > 12.5 among our z > 1.5 objects (low- 
est EWs) 1 _OursamDje also differs slightly from that 
of | Veilleux et al. (2009a) in that the strongest evolu- 
tion in the AGN fraction of local ULIRGs appears to 
be between phases III (early mergers) and IV (advanced 
mergers): the progression between phases IV and V (old 
mergers) being weaker. We observe, au contraire, only 
a small decrease in the PAH EWs of our phase IV ob- 
jects compared to those of phase III, and find most of 
the drop to occur after (rather than during) coalescence. 
We argue that the lack of a strong drop between phases 
IV and V in the 1-Jy sample is probably due to the se- 
lection at 60 jam which highly favors composite systems 
over AGN-dominated SEDs, and that this effect is absent 
from our sample because our selection is equally sensitive 
to both. Despite appearances, there is, thus, no clear in- 
dication that our objects differ much, intrinsically, from 
local ULIRGs in terms of their PAH properties along the 
merging sequence. 

Focusing in more details onto specific types of objects 
in our sample, we can draw the following conclusions. 
AGN-dominated spectra are prevalent in the later stages 
of the merging sequence, but can be found in almost 
all classes. Starburst-dominated spectra, on the other 
hand, disappear after the coalescence phase, and all reg- 
ular bulges are powered purely through AGN activity, 
consistent with the scenario in which star formation is 

ically carry much less flux than the 7.7/ira complex. In addition, 
the contribution of the continuum is also less important at high 
redshift due to bandpass shrinking. 
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Fig. 13. — Main panels: 11.3/xra and 7.7/xm equivalent widths of objects in our sample as a function of their redshift. Objects without the 
proper spectral coverage are omitted from the respective plots, but every object appears at least once. Insets: Distribution of equivalent 
widths split by redshift. The thick black histogram represents objects with z > 1.5, the thin grey one, objects with z < 1.5 and the dotted 
histogram, objects of unknown redshift. Objects with no detected PAH emission were all added to the lowest bin. Higher redshift objects 
tend to have lower PAH EWs in our sample. 
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Fig. 14. — Distribution of the 11.3/ira PAH equivalent width for z < 1.5 objects split into morphological classes. The dotted histogram 
represents the distribution of the sub-sample as a whole and has been slightly shifted for clarity. The grey area delineates the region of 
AGN-starburst composite systems. Starburst-dominated spectra fall to the right of that region whereas AGN-dominated ones populate the 
low-end of the distributions. Objects without a detected feature at 11.3/ira were all put into the lowest bin. The dashed lines represent 
the median EW of each morphological class. This figure demonstrates an average progression from star formation to black hole accretion 
along the merging sequence. 
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quenched in the process of bulge formation. Advanced 
mergers come in all three types of spectra, but are most 
highly peaked in the composite region, consistent with 
the idea that they are in a transition phase. 

Faint & Compact objects are interesting in that none 
of them have detected PAHs, yet they all clearly show 
low sersic index profiles, and whether their profile is due 
intrinsically to a rotating disk or to dust lanes obscuring 
part of a nucleus, both cases require significant quantities 
of gas. The reason why these objects do not show more 
star formation is therefore unclear. It is possible that 
feedback from the AGN acts to inhibit star formation. It 
is also possible, however, that star formation does occur, 
but that PAH features are either completely overpowered 
by the strong AGN continuum and/or that they are being 
destroyed by the radiation from the AGN. Alternatively, 
the star-forming environment in those galaxies could be 
so hot and dense as to wipe out any PAHs, or PAH fea- 
tures, itself, although that would be unprecedented. We 
show in the next section that these faint & compact ob- 
jects are actually highly obscured at mid-IR wavelengths, 
which implies a strong absorption of the radiation from a 
central hot source (AGN) by colder, obscuring material, 
or gas. This suggests that AGN radiation might well be 
responsible, perhaps in combination with a compact gas 
distribution, for modifying the PAH properties of these 
sources. 

4.6. Relation between optical depth and morphology 

The silicate absorption at 9.7 fim is a prominent sp ec- 
tral feature of many local ULIRGs ( |Hao et al.|[2QQ7| . It 
has been shown to arise whenever the main source of 
mid-infrared radia tion itself is veiled b y a column den- 
sity of colder dust ( |Levenson et al.|2QQ7| ) , and is therefore 
often interpreted as the signature oi an obscured AGN. 
The most buried of them are, in turn, seen as nascent 
quasars in the quas ar formation scenario proposed by 
Sanders et al. (1988). We compare, below, the morphol- 
ogy of objects with various degrees of obscuration. In 
particular, we examine how morphology can help us un- 
derstand the most obscured sources. We u se the rg^^m 
silicate optical depth, as defined in section |2.6.2[ to this 
effect. 

We first examine the relation of Tg^^ m with merging 
phase with the aid of Figure [T7| in which we illustrate the 
morphological composition ofobjects at different values 
of T9.7 Mm . Figure [l7| demonstrates that the highest op- 
tical depths are acETeved in the early phases of merging 
and through coalescence. Once objects reach late phases, 
their optical depth at 9.7/im diminishes considerably or 
vanishes altogether, as illustrated by the substantial in- 
crease in the fraction of late mergers at lower optical 
depths. Isolated objects also do not attain the large op- 
tical depths achieved in mergers. 

4.6.1. The Spoon Diagram 

More insight can be gained by looking at correlations 
with more than one parameter. One particularly useful 
diagram is that o f Tg^^ m versus PAH equivalent width. 
Spoon et al. (2007) showed that local ULIRGs lie on two 
distinct branches m that diagram: a horizontal branch 
of low optical depths spanning the whole range of EWs 
and a diagonal branch going from shallow silicate ab- 
sorption at high EWs to high values of Tg^^ m at low 
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Fig. 17. — Merger composition in four bins of Tg.7^ m . Only early 
and coalescing mergers populate the highest optical depths. Late 
mergers have predominantly low values of rg^m. 



EWs. They found the region comprised of objects with 
low EWs and medium optical depths, located in between 
the two branches, to be scarcely populated. These obser- 
vations led the authors to speculate that ULIRGs mainly 
evolve upward along the diagonal branch before quickly 
expelling their obscuring material to become quasars or 
Seyfert 1 galaxies on the low-EW end of the horizontal 
branch. Our morphological information allows us to test 
these ideas. 

We show the diagram of optical depth versus PAH 
equivalent width for low and high redshift objects in our 
sample separately in Figure [THJ Unfortunately, since our 
detectability is limited to PAH equivalent widths above 
0.1, we do not possess the sensitivity to probe into the 
AGN part of the diagram where the dichotomy between 
the two branches occurs in local ULIRGs. Nevertheless, 
we clearly detect a horizontal branch at low values of 
T9.7 Mm , so much that we show the morphological distri- 
bution of objects along that branch in binned histograms 
rather than as individual points, owing to their large den- 
sity. On the other hand, the number of objects with 
highly obscured AGNs (top-left corner of the diagram) 
is rather low in our sample, with only four such candi- 
dates (circled objects in Figure 18), insufficient to judge 



of the presence of a diagonal branch. 

Despite their small numbers, our morphological anal- 
ysis reveals that all of our buried AGN candidates have 
very specific host morphologies. Namely, they all belong 
to the Faint & Compact category. This observation indi- 
cates that the galaxies these objects live in possess large 
concentrations of gas in their center but do not have fully 
formed bulges. We speculate that these objects form in 
special merger configurations that result in the funneling 
of the major part of the gas and dust content into the 
center without much loss to stripping and tidal features. 
We think that such mergers might occur in low speed 
encounters. 

Furthermore, the highly obscured nature of these ob- 
jects implies that they have not been able to clear the gas 
surrounding their central hot source, or active nucleus, 
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Fig. 18. — Spoon diagram jSpoon et al. 2007) for our sample color-coded by morphology: cyan for isolated spirals, blue for early mergers 
and close pairs, red for advanced mergers, green for faint & compact objects, and yellow for late mergers. Owing to the large density of 
points at Tg^^m < 1.12, we show histogram distributions in bins of EW, rather than individual points. Objects with upper limits on their 
EW were all placed in the AGN bin. Objects of unknown redshift were split half-and-half between the two diagrams. They all fall in 
the lower-left bin, by definition. Bars of the histogram run fro m to 10 objects. Circled objects represent obscured quasar candidates. 
Their morphology, IRS spectrum and SED are shown in Figure [l9| They all possess faint & compact morphologies. Late mergers tend to 
aggregate in the lower-left part of the diagram where unobscured AGNs live. Other morphological classes do not occupy preferentially any 
specific region. 



yet and would therefore have to be very young mergers, 
i.e. still in the early coalescence phase. We would expect 
the stellar component to begin to be randomized, and 
the smooth exponential profile to start to be destroyed, 
at that stage. The observed sersic indices (n < 2) of these 
objects are consistent with such a process beginning to 
take place. These objects thus appear to be growing 
their bulges and black holes simultaneously and to be 
perfect examples of bulge/black-hole co-e volution and of 



the c lassical quasar formation scenario of |Sanders et al. 
([1988) 



Most mergers, however, do not go through a faint & 
compact phase. In fact, most objects going through coa- 
lescence show tidal tails and/or streams, as well as higher 
sersic indices, and thus belong to the advanced merger 
class. Having already built most of their bulge, these 
objects will never evolve into faint & compact objects. 
Because of their tidal features, these objects are also un- 
likely to have evolved from a faint & compact phase, thus 
making the two evolutionary paths largely independent. 
Advanced mergers show, on average, shallower silicate 
absorption, as illustrated in Figure 18 Rather than go- 
ing through that obscured quasar phase, most infrared- 
luminous mergers thus appear to either move off and on, 
or about, the horizontal branch instead, with often some 
combination of both PAHs and silicate absorption. 

Conversely, we know of at least one example of a high 
7"9.7/im and weak PAH source locally that is not a faint 
& compact source. This is IRAS08572+3915, an early 
merger and one of the objects we simulated high redshift 
observations for. Taken individually, however, both in- 
teracting galaxies in this source loose their tidal tails and 
extended features at high-z and appear faint & compact. 
Moreover, the two bodies are only 6 kpc apart and are 
likely on the brink of merging and becoming a true faint 
& compact source. Nevertheless, this example illustrates 



that the obscured quasar phase can sometimes occur be- 
fore final coalescence. 

As for early and advanced mergers, Figure [18] shows 
their distribution in the Spoon diagram to be fairly well 
mixed implying little evolution between the two phases. 
The large parameter space spanned by both those cate- 
gories further suggests that there are many possible paths 
a merger event can traverse. It appears that the only pre- 
dictable stage is the end-point of the merging sequence: 
the late mergers, shown to aggregate in the lower left 
corner of the diagram, where both star formation and 
obscuration are low. 

Our sample, thus, does not attest to any precise evo- 
lutionary path in Tg m 7 firn , - EW (PAH) space, such as sug- 
gested by Spoon et al. ( |2QQ7| . Rather, it seems that 



many avenues are possible. Nevertheless, it appears that 
most objects gravitate around the horizontal branch, and 
only a certain kind of mergers make it to the top-left 
corner: those with very dense and dusty cores that are 
classified as faint & compact at high redshift. However, 
all objects appear to end up in the low-rg.7 Mm , low-EW 
corner at the end of the merging process. 

5. DISCUSSION 

We have shown that a high fraction, about 80%, of 
bright {S{2Ajim) > 0.9 mJy), high-redshift 24/im-sources 
are likely to be ongoing mergers in various stages of the 
process. We further demonstrated that although our 
sources show a broad range of morphologies and populate 
the entire merging sequence, the relative number of early 
mergers increases substantially at z > 1.5. We showed 
that star formation activity as probed by the equivalent 
width of PAH features, decreases, on average, along the 
merging sequence to the profit of black-hole accretion. 
We finally demonstrated that obscured quasars, in our 
sample, live in "faint & compact" galaxies. 
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Fig. 19. — NICMOS images, residuals after subtraction of the smooth profile, and IRS spectra (smoothed in red) of all obscured 
(T9. 7^771 > 3.36) quasar candidates in our sample. Images are shown in negative, whereas residuals in positive. Dotted lines in the spectra 
show the location of PAH features. The arrow marks the center of the 9.7/ira silicate absorption feature. All objects are faint &; compact 
and three out of four show faint residuals, the shape of which is characteristic of mergers in two of the cases. 



We mirrored our analysis on that of Veill eux et al. 
fl2002| |2QQ6[ |2QQ9a| on the 1-Jy sample, and simulated 
high-redshift JNICMOS observations of local ULIRGs in 
order to address the question as to whether our sources 
are high-redshift analogs of local ULIRGs. To a large 
extent, we presented our comparison of the two samples 
concurrently with our results in section 4, but discuss 
that question in a broader context in this section. The 
1-Jy sample, however, is arguably not the best point 
of comparison, since the sample is selected at a sub- 
stantially different rest-frame wavelength. Furthermore, 
comparison with local galaxies gives us only part of the 
picture; comparison with other high-redshift samples is 
also highly desirable. We propose to carry such compar- 
isons in this section. We will discuss our results in the 
context of five types of infrared-luminous objects: lo- 
cal ULIRGs (from the 1-Jy sample), other bright 24/im- 
selected samples at high-redshift, sub- mm galaxies, and 
finally z ~ 1, 70/im-selected objects. Lastly, we com- 



pare our results with theoretical predictions and discuss 
implications for our understanding of galaxy evolution. 

5.1. Are bright high-redshift 24/ira galaxies analogs of 
local ULIRGs? 

Our observations suggest that, up to z = 1.5, 
they are analogs of, not only local ULIRGs, but 
also top-end LIRGs (L IR > 1O 115 L ). This is evi- 
denced by their large infrared luminosities (£3-1000/1™ ~ 
10 11,5 -10 12,5 L©) and their high merger fraction (up to 
83% at z < 1.5), which are the two defining c haracter- 
istics of high-end LIRGs and ULIRGs loca lly ( |Murphy 
et al.|1996||Veilleux et al |2QQ2|[lslnda|2QQ4[ ). Our analy- 
sis, however, has revealed that our 24/ira-selected sample 
at z < 1.5 tends to draw from objects in more advanced 
stages of merging, likely owing to the shorter wavelength 
selection. When compared to redshifted galaxies from 
the 1-Jy sample, our objects also tend to show stronger, 
brighter merger signatures such as more visible tidal fea- 
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tures and a higher incidence of detected residuals. Spiral 
galaxies make for only 17% of our objects at z < 1.5, 
thus rejecting the hypothesis that our sources are scaled 
up versions of lower luminosity local LIRGs (which hav e 
disk morphologies in > 40% of cases; [Wang et aL| |2006), 
even though some of our spirals do show slightly elevated 
IR-luminosities (up to 1.6 x 10 12 L©). 

At redshifts above 1.5, we observe a dramatic rise in the 
number of early mergers among our bright 24/im sources 
despite the fact that we are probing, at those redshift, 
only extremely high luminosities {Ljr > 1O 12,5 L ), and 
despite the fact that > 90% of these objects have PAH 
EWs indicative of the presence of an AGN (EWr^m < 
1.2). Locally, this combination of high luminosity and 
low EWs can only be found among more advanced merg- 
ers. These early mergers that represent, at z > 1.5, half 
of our sample, thus form a new type of sources. Morpho- 
logically, albeit very messy, they do not appear funda- 
mentally different from lower redshift luminous mergers 
at similar stages. Their PAH EWs, however, indicate 
that they are primarily AGN-dominated or composite 
systems. It is, therefore, an earlier triggering, and sus- 
tained fueling, of a strong obscured AGN that distin- 
guishes them from local ULIRGs. As we discuss in sec- 
tion |5.5[ we think this might be the consequence of the 
elevated gas content o f these objects in comparison to 
that of local ULIRGs flYan et al".||20To| >. The infrared- 
luminous aspect of galaxy evolution thus appears to pro- 
ceed differently at high redshifts, in a way that is not 
well represented by local ULIRGs. 

5.2. Comparison with other samples of 24/im -selected 

galaxies 

Dasyra et al. (|2008[) presented a first analysis of a 



z > 1.5 subset ol the sample used in this paper that 
was selected from objects in our first Spitzer /IRS pro- 
gram. By looking for merging pairs, they were able to 
place a lower limit on the number of interactions in z ~ 2 
bright 24/im galaxies at 52%. We have expanded those 
results by performing a systematic search for merging 
signatures and find that the other half of our z > 1.5 
objects all show signs of possibly being later-stage merg- 
ers (from coalescence onward). 60% of those objects are 
identified with high confidence, being of confidence level 
4 or less (cf. § 2 5.2). 

|Dasyra et al.| ( |2008| ) further found hints that these 
high-redshilt Brignt 24/im sources might be disk- 
dominated rather than bulge-dominated like local 
ULIRGs. Our analysis reveals that this is indeed the 
case and that two thirds of our objects at z > 1.5, 
in fact, possess disk-dominated profiles (B/D < 0.5 or 
n < 2.3). This can be largely attributed to the high frac- 
tion of early mergers at those redshifts, but many, more 
advanced, mergers also show rather low B/D ratios or 
sersic indices. We speculate that this is due to a combi- 
nation of two effects: high gas fractions in those objects 
that act to maintain a disk around the central o bject for 
longer (and perh aps indefi nite) periods of time (|Robert- 



1000 colors. Bussmann et al. (2009) presented NIC- 
MOS (as well as ACS/WFPC^) morpEologie s of 31 such 
sourc es selected from the larger sample of |Dey et ah] 
(2008 ) and chosen to have 24 /im fluxes in excess of 0.8 
mJy. [Melbourne et al.| ( |2009| ) presented NIR morpholo- 
gies of 15 additional DOGs obtained through if-band AO 
observations, eleven of which have 24/im fluxes above 0.8 
mJy. Both papers find DOGs to possess a high fraction 
of disk-like profiles (~ 90% and 50%, respectively), and 
argue, based on their size, sersic index, axis ratio and 
morphological parameters that DOGs are consistent with 
being late- stage mergers, transitioning from the chaotic 
coalescence towards becoming relaxed ellipticals. 

There are 35 DOGs in our sample (Figure [T]), 75% of 
which are at z > 1.5. Not counting z < 1.5 objects, we 
observe that about half of our DOGs are actually pre- 
mergers (phase III objects), while only half are in more 
advanced stages, contrary to the hypothesis that they 
are predominantly late-stage objects. We agree, however, 
that most of them have disk-dominated profiles, and find 
that fraction to be 2/3. The morphological properties of 
DOGs, summarized in Table [3[ suggest that they are not 
very different from the overall 24/im-bright population 
at z > 1.5. They only tend to draw a little bit more 
towards AGN-dominated morphological classes, that is 
faint & compact objects and late-stage mergers. As a 
consequence, we also find them to be slightly more com- 
pact, AGN-dominated and obscured than average bright 
24/im- selected galaxies at z > 1.5. 

5.3. Comparison with sub-mm selected galaxies 

Sub-millimeter selected galaxies (SMG s) form another 
type of galaxies known to lie at z « 2 ( |Chapman et al. 
2005[) and to be ex tremely luminous in the infrared 
( |Kovacs et al. 2006). These objects are thus closely 
related to our bright 24/im-sources. Unlike our 24/im- 
sources, however, SMGs po ssess PAH- dominated mid-IR 



spectra in > 80% of ca ses ( [Pope et al.||2008 Menendez- 
|Delmestre et al.| |2009[ ) . They would, therefore, largely 
populate the high-EW, high-z part of Figure 13j the 
only region not occupied by our galaxies. The two pop- 
ulations, thus, complement one another in covering the 
broader high-redshift ULIRG population. Nevertheless, 
many SMGs show composite spectra, so that there is 
also significant overlap between the two types of sources 
(Sajina et al. in preparation). 

High-resolution HST observat ions of sub-millimeter 
galaxies (SMGs) a t both optical |Sniail et al.|1998 



selice et al.l|2003l jPope et al.|[200b| |L)unlop eFaT 
Swinbank et al.| |201 0D and near- infrared wavelengths 



Con- 



( Swinbank et al. 2010) have revealed that a large fraction 
ol them (407o-907o) show highly disturbed or multiple 
component morphologies indicative of merging activity. 



ture of a majority of these objects ( 


Neri et al.|2003| Greve 


et al. 


20051 Tacconi et al. |2006, 5 


?008; Schinnerer et al. 


2008| 


lono et al. 2009| Bothwell et al.||2010p, although 



large star- forming disks have certainly also been found 



|son et aL||2006| see also § 5.5), and band-shifting that to exist among the sub-millime ter population (Carilli 



causes us to see those objects m their rest-frame R, V or 
B band where their disks can be readily detected. 

Another type of bright 24/im-selected galaxies are the 
so called DOGs (for dust-obscured galaxies). These are 
z ~ 2 objects with very red observed F(24fj J m)/F(R) > 



et al.||2010| |Bothwell et~aL 



vor a me rger origi n for Slv 



2010]). Most models also la- 
Gs 



(e.g. 



jBaugh et al.||auub 
Chakrabarti et aL||2QQ8| |Narayanan et al. |2UU9h2(JlUbp 
even though regular disk star formation fueled b y cold 
flow accretion also appears to be a viable path (Dave 
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TABLE 3 

DOG MORPHOLOGIES COMPARED TO 24/im-GALAXIES a 



Merging sequence Morphological class 



All z > 1.5 
24/ira-galaxies 



DOGs 



Isolated objects 



All spirals 



First approach (Jlose pairs (phase 1) 



Fariy mergers 



First contact (phase 11) 'Z%(l)~ 
rre-m ergers (phase ill) 

Triplets 67o ( 'I ) ' 



U 



4t>7o ) Wo [12) 



Advanced 



Uoalescence 



Mergers (phase IV) 
Fain t &; Com pacE 



2(6% ( lb. 4') 
V% ('4.4) ~ 



iy% tb 1 ) 



V2% C3) 



Fate mergers 



Uld mergers (phase V ) 
rtegular bulges 




a After red istribution of objects with unknown redshifts in the proportions 
found lor objects with known redshifts. 



et1^|2QlQ ). 

For our purpo s es, we mostly refer to the results of 
Swinbank et al. (2010), who published the first mor- 
phological analysis ol sub-millimeter galaxies at near- 
infrared wavelengths. Although they did not use the 
same approach towards morphological classification as 
we do here, we can still infer, from their i^-band im- 
ages, that 40% of objects in their sample have two or 
more distinguishable components. This makes the frac- 
tion of early mergers among SMGs very close to that of 
our z > 1.5 galaxies. We can also see that their objects 
display many of the same kind of disturbed morphologies 
as ours. We might, thus, expect the morphological dis- 
tribution of SMGs not to be hugely different from that 
of our high-z objects, especially given that CO studies, 
as mentioned above, have found a number of those ob- 
jects to possess rather chaotic kinematics characteristic 
of the more advanced merging stages. Nonetheless, we 
know that purely AGN-dominated objects are not found 
in sub- millimeter samples, and, based on the results pre- 
sented in section |4.5[ we would therefore expect faint & 
compact objects, as well as very late-stage objects that 
appear as regular ellipticals, to be absent from SMG sam- 
ples. This does appear to be true though more detailed 
analysis would be needed to confirm it. Size and struc- 
tural analysis, which we postpone for subsequent publi- 
cation, might also reveal more subtle differences between 
24/im and sub- millimeter galaxies. 

5.4. Comparison with far-IR selected galaxies 

Kartaltepe et al. (|2010|) have conducted the most com- 



prehensive study of the role of mergers in far-IR selected 
galaxies to date using the full 2 sq. deg. of HST/ACS 
imaging and multi-wavelength coverage in the COSMOS 
field to investigate the morphological properties of a com- 
plete sample of 1500 70/im-selected galax ies. They too 
adopted the morphological classification of |Surace] (|1998) 
based on the merging sequence, enabling direct compari- 
son with our results. This comparison, however, is a little 
bit complicated by the fact that they use optical rather 
than NIR imaging, causing ~ 20% of their objects at 
z > 1 and Ljr > 1O 12 L to be too faint to classify mor- 
phologically. This caveat is sufficiently minor though to 



be circumvented, and we demonstrate below that inter- 
esting conclusions can still be drawn from a comparison 
of the two samples. 



Visually, the 70/im-galaxies of Kartaltepe et al. (2010) 
resemble a lot more our low-redshilt (z < 1.5) sample 
than either our high-redshift galaxies or the SMGs. Since 
their sample includes many low-redshift objects with 
moderate IR-luminosities, we use, for our comparison, 
only those objects in their sample with a log(L/^/L ) > 
11.5. We also limit our sample to z < 1.5 objects in 
order to ensure that both samples have comparable dis- 
tributions in both redshift and total IR-luminosity. Ta- 
ble [4] shows the detailed com parison by morphological 
class. In their classification, Kartaltepe et al. (2010) 
make the distinction between minor and major mergers 
whereas we do not. Then, many of their sources have 
unknown morphologies because they are too faint to be 
detected in the observed optical. For comparison pur- 
poses, we, therefore, redistributed their minor mergers, 
for lack of a better prior, between phases III & IV in the 
same proportion as that of their major mergers, and their 
unknown objects equally among phase III and "faint & 
compact" objects, since these are the two classes our faint 
objects fall in. We show in table [4]both the observed and 
redistributed numbers. 

Table [4] confirms the astonishing correspondence be- 
tween our z < 1.5 sample and their bright 70/im-sources. 
It is remarkable that the proportions in all phases of the 
merging process differ by no more that ~ 15% (of the 
total sample), lending great support to both our results. 
Nevertheless, the increase in the number of isolated spi- 
rals and early-phase mergers among 70 fim sources and 
the corresponding decline in the number of coalescence 
and post-coalescence mergers when compared to 24/im 
sources is clear, and probably a simple consequence of the 
higher proportion of star-formation dominated galaxies 
among 70/im-selected sources. 

The broad correspondence between 24/im and 70/im- 
selected galaxies suggests, as does the similarity between 
z ~ 2, 24/im- selected and submillimeter-selected galax- 
ies, that most of the observed morphological specificity 
in these different samples is common and intrinsic to all 
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TABLE 4 

Comparison of morphologies of 24 and 70/xra-SELECTED galaxies 



Merging sequence 


Morphological class 


24/ira-selected a 


70/ira-selected b 
(observed) 


70/ira-selected c 
(estimated) 


Isolated objects 


All spirals 


16% 


22% 


22% 


First approach 


Close pairs (phase I) 


4% 










First contact (phase II) 


3% 


1% 


1% 


Early mergers 


T~» „ „ „ / 1 TTT\ 

Pre-mergers (phase III) 


i c 07 
lO70 


01 07 

21% 


36% 




and Triplets 










Advanced 








Coalescence 


Mergers (phase IV) 


34% 


15% 


22% 




Faint & Compact 


9% 




5% 




Old mergers (phase V) 


15% 


3% 


3% 


Late mergers 


Regular bulges 


4% 


9% 


9% 




Pure point sources 


1% 


2% 


2% 




Minor mergers 




16% 






Unknowns 




11% 





a z < 1.5 obj 
b Data from 
c After redis' 



ects only 



Kartaltepe et al. (2010), including only objects with log(Lm/ Lq) > 11.5 
ribution ot minor mergers and unknowns. 



infrared-luminous objects, largely irrespective of whether 
they are powered by star formation or an AGN (although 
these can be responsible for small differences between the 
various classes). It also suggests that this commonality 
springs from the origin of those objects in massive merg- 
ers. The difference between z ~ 1 and z ~ 2 samples, on 
the other hand, indicates that there is redshift evolution 
in the morphological character of IR-luminous objects. 
We argue in the next section that this difference is the 
result of higher gas fractions at higher redshifts. 

5.5. Comparison with theory and consequences for 
galaxy evolution 

Hydrodynamical simulations of merging disk galaxies 
combined with radiative transfer calculations have pro- 
vided us with a succe ssful physical model of the origin of 
ULIRGs an d QSOs flDi Matteo et al.||2QQ5| |Narayanan 



et al.|2010a ) . These can be used, more generally, to study 
the outcome of any merger by varying the set of initial 
conditions given by parameters such as the mass ratio 
of the two galaxies, their gas fractions, their relative ori- 
entation and so on. Using the statistical results of such 
a suite of simulations, and embedding them into a cos- 
mological framework using a n observational ly motivated 
halo occupation distribution, |Hopkins et al.[ in a series of 
papers starting from (2008), have been able to calculate 
the role mergers play in the evolution of various global 
and observable p roperties of the galaxy population. 

Most recently, [Hopkins et al.| ( |2010[ ) calculated the 
IR-luminosity at which, according to their models, star 
formation transitions from occurring mostly in isolated 
galaxies to being primarily triggered by mergers. We 
plotted their results on top of our Figure [IT] for compar- 
ison. Our data agree with their models on that we find 
nearly all of our isolated spirals below the predicted line, 
indicating that their models provide an accurate upper 
limit as to the luminosities at which quiescent star forma- 
tion usually occurs. We also find good agreement at high 
redshift where nearly all of our objects fall above their 



line and are, indeed, found to be mergers. We appear 
to disagree, on the other hand, at lower redshift where 
we find many merger-identified objects below their di- 
viding line. This, however, is not necessarily a problem, 
the reason being that the 24/nm selection of our sam- 
ple is very different from that of an Ljr selection. In 
particular, we know that our selection highly favors ob- 
j ects whose SEP contains an obscured AGN component 
( |Sajina e t al. 2007; Dasyr a et al ■ 112009 D . If we adopt the 
results ot Hopkins &; Quataert (2010) which argue that 
the fueling of a strong AGN (i.e. M > 0.1 M yr" 1 or 
Lir > 10 11 L yr _1 ) can only occur in mergers, then 
this means that we are preferentially selecting mergers, 
and that our sample is biased compared to the general 
infrared-luminous population. Our data, thus, do not 
lend th emselves very well to a comparison with the re- 
sults of |Hopkins et al. (2010) except on isolated objects, 
as discussed above. They do, however, seem to support 
another result, n amely the one argued for in Hopk ins fc| 



Quataert 
strong AG 



(2010) which is that it is very hard to fuel a 
JIN other than through a merger event. 
Another aspect of those galaxy merger simulations is 
that the merging of two massive gas-rich disks almost 
always produces an extended IR-luminous phase domi- 
nated, at first, by star formation during early stages of 
the merging, and then by accretion onto the central black 
hole after coalescence. Our res ults, as well as those from 
observations of local ULIRGs ( Veilleux et al.|2009b 



Far- 



|rah et aL]|2009| ), show that the above scenario is, indeed 
likely to be the evolutionary path of many low to medium 
redshift (z < 1.5) ULIRGs. However, both our observa- 
tions and that of local ULIRGs also show that there is a 
great variety of objects at all stages of the merging pro- 
cess, implying many more possible evolutionary paths. 
This variety has not been fairly reproduced, so far, in 
simulations. 

At high redshifts (and high luminosities), our results 
indicate that mergers depart even more from that classi- 
cal scenario, and that black hole accretion plays an im- 
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portant role, at least among our sources, already in the 
early stages of merging. CO observations have demon- 
strated that 24/im-selected galaxies at z ~ 2 possess 



larger amounts of gas than do local ULIRGs ( |Yan et al. 
2010). We speculate that these elevated gas tractions 
might be at the origin of the increased AGN activity at 
early stages. 

Another characteristic of bright 24/ira-galaxies at z ~ 
2 that our observations, as well as that of others (e.g. 
|Dasyra et al.|2008} |Bussmann et al.|2009 ), have revealed 
is their preponderantly disk-dom inated profiles. Here, we 
know, thanks to simulations by Robertson et al.| ( |2006[ ) 
and others, that high gas fractions m major mergers can 
act to stabilize the disk component and that major merg- 
ers can even result in rotationally supported disk rem- 
nants provided that the prog enitors are extremely gas- 
rich (f q0iS > 0.8) (see also ISpringel & Hernquist||2005[ 



Barnes||2002| ). The results oi these simulations thus pro- 



vide a natural explanation linking this large observed 
fraction of disk-like morphologies to the high gas content 
o f these objects. 
Robertson et al. (2006|) argue that the kind of galaxy 



assembly in which mergers result in a disk remnant is 
probably limite d to high redshifts wher e gas reservoirs 
are much larger ( |Noterdaeme et al.|2009 ). Although it is 
impossible to know without detailed kinematics whether 
any of our z ~ 2 mergers exhibit such behavior, we 
see, based on their profile, visual appearance and merger 
stage, two such candidates: MIPS16144 and MIPS16122. 
Both have clearly detected PAHs (EW(11.3/xra) = 0.5 
and 0.19, respectively). 

These two tentative disk merger remnants at z ~ 2, 
and the lack thereof at lower redshifts, combined with the 
overall large fraction of objects with disk profiles at those 
redshifts, suggest that this epoch might correspond to the 
end of the era of gas-rich mergers a la |Robertson et aT7| 
( |2006| ). Interestingly, that period in cosmos history also 
concurs with the onset of the rapid r ise in mass density of 
elliptical and red sequence galaxies ( |Arnouts et al.| 2007 



Ilbert et al.|2010[ ). These two observations can naturally 
explain each other, and point to a transition, at z ~ 2, 
from an epoch of disk-dominated galaxy evolution to an 
epoch of bulge formation. 

6. SUMMARY 

Our new analysis of 134 HST/NICMOS images of 
bright high-redshift 24/im-selected galaxies reveals a 
higher incidence of mergers than previously reported. 
Our study is based on a flux-limited sample of objects 
selected at 5(24/ira) > 0.9 mJy, with a sampling of the 
[24] — R color space that favors high redshift objects. 
Our sample shows two broad peaks in its redshift distri- 
bution, one at z ~ 1 and one at z ~ 2. Full SED analysis 
presented in a companion paper (Sajina et al., in prepara- 
tion) demonstrates that most of our objects above z = 1 
have total IR-lumonisities above 1O 12 L , whereas those 
at lower redshifts mostly have luminosities ranging from 
10 11 ' 5 to 1O 12 L . Above z = 1.5, nearly all objects in our 
sample have infrared luminosities in excess of 1O 12,5 L . 
The merging properties of our sample are summarized 
below: 

• We find a high overall merger fraction, which we 
estimate to lie at ~ 80%. Possible values, however, 



range from 62 to 91%. 

• Mergers represent 65 to 96% of ULIRGs in our sam- 
ple, but our best estimate puts that fraction at 
87%, suggesting little evolution in the dynamical 
origin of that category of objects. 

• Isolated spirals form 9% of our sample. They all 
lie at z < 1.2 and have luminosities of £3 1000/2™ < 
1.6 x 1O 12 L . 

• At z < 1.5, the distribution of our objects draws 
towards more advanced stages of the merging pro- 
cess when compared to local ULIRGs. Both sam- 
ples, however, are dominated by singly-nucleated 
(coalesced/coalescing) mergers with bulge-like pro- 
files. 

• Comparison with simulated observations of red- 
shifted local ULIRGs demonstrate that our ob- 
jects also tend to show stronger merger signatures. 
This manifests itself in brighter tidal tails and a 
higher detection of residuals after subtraction of 
the smooth component. 

• At z > 1.5, the morphological composition changes 
significantly with a dramatic rise in the number of 
early mergers: systems in which the two merging 
objects are still distinguishable. 

• Still at z > 1.5, 60% of our objects possess disk- 
dominated profiles which we argue is partially a 
consequence of their higher gas content that acts 
to retain a larger disk for longer periods of time 
during a merger. 



Mid-IR I RS spectra for all ob ject s in our sample wer e 
presented in Sajina et al. Q2007 ) and Dasyra et al. (2009). 
These data give us an indication of the power source at 
the origin of the bright mid-IR flux of our objects in the 
form of the equivalent width of the PAH features they 
probe. They also provide a measure of AGN obscuration 
through the depth of the 9.7 [im silicate feature. Combin- 
ing this information with our morphological data gives us 
new insights on the merging process and the formation 
of quasars at high redshift. It allows us to deduce the 
following conclusions: 



• Statistically, the importance of AGN activity rela- 
tive to star formation increases as objects progress 
along the merging sequence. At z < 1.5, most ob- 
jects are dominated by star formation before their 
first pass, whereas past coalescence, most become 
AGN-dominated. This t rend is similar to that ob- 
served in local ULIRGs flVeilleux et~aT|2009a| ). 



• At z > 1.5, however, even early mergers possess, in 
> 90% of cases, a significant power-law continuum 
indicative of the presence of an AGN. Combined 
with the fact that these have luminosities in ex- 
cess of 1O 12 5 L , and the rarity (or inexistence) of 
such objects locally, this signifies that black- hole 
accretion plays an important role earlier on in the 
merging process at higher redshifts. We speculate 
that this, again, might be a consequence of the el- 
evated gas fractions of objects at those redshifts. 
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• We find that obscured quasars live in hosts with 
faint & compact morphologies. 



• Our simulated observations indicate that local 
dense-core mergers that have diffuse tails and not 
yet fully formed bulges are the best candidate 
analogs of faint & compact objects at high red- 
shift. These kind of mergers probably occur only 
in certain types of initial configurations. 

• We find, however, that only few objects, four in to- 
tal in our sample, possess spectra characteristic of 
obscured quasars. We conclude that either merg- 
ers rarely go through that phase or else it's very 
short-lived. 



well those of 70/im-selected galaxies of similar luminos- 
ity and redshift. At z > 1.5, we begin to see a new 
type of galaxy appearing in large numbers: early merg- 
ers with high infrared luminosities (Lm > lO 12,5 I/ ) and 
low PAH equivalent widths (EWV^m, ~ 0.3). These 
have no known local equivalent. We also begin to see 
a preponderance of disk-dom inated profiles, which con- 
trasts that of local ULIRGs flVeilleux et"aL]|2006| >. On 
the other hand, these high-z bright 24/ira-galaxies ap- 
pear not to differ as much from other types of z ~ 2 



infrared-lu minous objects such as DOGs ( |Dey et al. 2008) 
or SMGs flBlain et aL]|2002| ). We therefore conclude 
that, although still driven by mergers, high redshift ultra- 
luminous infrared galaxies are substantially different in 
nature from local ULIRGs and that this difference i s 
likely rooted in their higher gas masses ( Yan et al.|2010). 



• We find late-stage mergers to predominantly show 
spectra characteristic of unobscured AGNs. 

• We do not find other morphological classes to oc- 
cupy any definite region in the Spoon diagram, im- 
plying that infrared-luminous mergers do not fol- 
low any characteristic trajectory, but can evolve 
through many possible paths. 

Comparison with other galaxy samples indicates that 
up to z ~ 1.5 the properties of bright 24/im-selected 
galaxies are very similar to that of local ULIRGs and 
top-end LIRGs (L IR > 1O 115 L ). They also match 
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APPENDIX 
REDSHIFTING PROCEDURE 

We describe, below, the procedure we followed to create simulated observations of local ULIRGs and mergers at 
high redshift. The simulation procedure involves two main steps. The first one consists of bringing the image to the 
NICMOS resolution at the desired redshift. To achieve this, we effectively resample the WFPC2 PSF on a grid whose 
pixel size equals ^'xg^ 1 ^ times the original pixel size of the PSF («s NIC2 corresponds to the pixel scale of the NIC2 
camera in arc seconds, or 0.0756" , and <s WFPC2 to that of the appropriate WFPC2 camera in arc seconds, or 0.1" for 
the WF cameras, and 0.0455" for the PC camera). This gives us the intrinsic resolution of our object at that redshift. 
We then deconvolve the NICMOS PSF with that resampled PSF by simply dividing in Fourier-space and applying a 
tapering at the edges, if necessary. This step is the one that involves the most user interaction in order to make sure 
that the division in Fourier-space does not create too much ringing, or does not overly amplify the noise. 

Once we have our convolution kernel in hand, we proceed to resample our image on a grid whose pixel size is 
multiplied by the same factor as above, and then convolve it with the WFPC2 — » NIC2 kernel. We then dim our 
image by the square of the ratio of the appropriate luminosity distances, plus a factor of (1 + z) to account for the fact 
that we are working with fluxes per unit wavelength, rather than bolometric quantities. We then further multiply our 
image by the ratio of the sensitivities of the two cameras, expressed in erg cm -2 A -1 / DN, and the ratio of exposure 
times, in order to get an image in units of NICMOS counts per second. We finally put that image onto an actual, but 
blank, NICMOS image obtained during the course of our observations. 

We repeated that redshift ing procedure four times, shifting the object around using the same dither pattern we used 
for our observations. We also added a random jitter component to our shifts to reproduce observations most faithfully 
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(though it turned out not to make much difference). Each image was laid onto a different blank background image. We 
then drizzle the four simulated images back onto one grid using the same procedure as for our actual data to produce 
our final simulated images. 

MORPHOLOGICAL DATA 



TABLE 5 
Source morphology 



MIPS ID 


Redshift 


m H 


t-i /o fknc^ 
'1/2 V^FW 


Dominant 


Morphological 


Confidence 








component 






MIPS34 


0.65 


17.84 


2.55 


bulge 


Old merger (V) 


4 


MIPS42 a 


1.95 






disk 


Faint & Compact 


6 


MIPS78 


2.65 


21.96 


2.84 


disk 


Faint & Compact 


4 


MIPS159 




18.43 




ambiguous 


Point source 


6 


MIPS168 


0.24 


17.61 


1.20 


bulge 


Advanced merger (IV) 


2 


MIPS180 


2.47 


22.08 


2.08 


disk 


Pre-merger (III) 


4 


MIPS184 




18.12 


~ 5.81 b 


disk 


Edge-on spiral 




MIPS213 


1.22 


20.61 


1.37 


disk 


Pre-merger (III) 


4 


MIPS224 


1.47 


18.55 


2.83 


bulge 


Advanced merger (IV) 


2 


MIPS227 


1.63 


18.65 


2.64 


bulge 


Advanced merger (IV) 


3 


MIPS268 


1.69 


19.90 


1.95 


disk 


Pre-merger (III) 


4 


MIPS289 


1.86 


18.27 


3.54 


disk 


Pre-merger (III) 


1 


MIPS298 


3.49 


22.08 




disk 


Pre-merger (III) 


5 


MIPS322 




19.38 


~ 0.96 b 


( t 1 1 1 l M ^ l l V M i 


Pre— merger (III) 


2 


MIPS324 


0.95 


18.40 


3.36 




Advanced merger (IV) 


3 


MIPS331 


1.03 


18.66 


0.90 


ambiguous 


Old merger (V) 


5 


MIPS351 


1.16 


19.44 


1.41 


^disk* 3 


Old merger (V) 


5 


MIPS358 


0.81 


18.75 


3.09 




Advanced merger (IV) 


2 


MIPS369 


0.70 


19.05 


3.28 


disk 


Face-on spiral 




MIPS397 


1.35 


19.32 


3.47 


bulge 


Regular bulge 


6 


MIPS419 


0.83 


19.61 


1.51 


disk 


Faint & Compact 


5 


MIPS446 


0.82 


18.37 


3.84 


ambiguous 


Advanced merger (IV) 


3 


MIPS464 


1.85 


20.86 


1.92 


disk 


Close pair (I) 


4 


MIPS488 


0.69 


18.13 


3.70 


disk 


Advanced merger (IV) 


2 


MIPS495 


0.75 


19.73 


3.00 


ambiguous 


Old merger (V) 


5 


MIPS505 


1.65 


21.12 


0.85 


bulge 


Advanced merger (IV) 


4 


MIPS506 


2.47 


20.78 


1.71 


disk 


Pre-merger (III) 


4 


MIPS512 


0.99 


19.33 


1.98 


bulge 


Advanced merger (IV) 


4 


MIPS544 


0.98 


19.68 


1.76 


disk 


Faint & Compact 


5 


MIPS546 


1.07 


18.55 


4.01 


disk 


Advanced merger (IV) 


3 


MIPS562 


0.55 


17.85 


6.88 


disk 


First contact (II) 


2 


MIPS7985 


2.78 


18.42 


1.26 


disk 


Advanced merger (IV) 


3 


MIPS8069 


0.70 


17.50 


1.68 


bulge 


Advanced merger (IV) 


4 


MIPS8071 


0.98 


18.42 


0.53 


bulge 


Regular bulge 


6 


MIPS8098 


1.07 


18.46 


3.19 


bulge 


Advanced merger (IV) 


4 


MIPS8107 


0.94 


19.02 


5.64 


disk 


Advanced merger (IV) 


3 


MIPS8172 


1.11 


19.10 


0.41 


ambiguous 


Old merger (V) 


5 


MIPS8179 


0.59 


18.04 


3.17 


bulge 


Advanced merger (IV) 


3 


MIPS8185 




19.29 


~ 0.58 b 


bulge 


Old merger (V) 


5 


MIPS8204 


0.85 


18.83 


1.75 


disk 


Pre-merger (III) 


4 


MIPS8224 




19.89 


~ 1.31 b 


bulge 


Old merger (V) 


4 


MIPS8226 


2.10 


19.86 


0.17 


bulge 


Advanced merger (IV) 


4 


MIPS8233 


0.99 


19.09 


1.92 


disk 


Pre-merger (III) 


4 


MIPS8242 


2.45 


19.72 


4.92 


disk 


Pre-merger (III) 


3 


MIPS8245 


2.70 


22.37 


1.63 


disk 


Faint & Compact 


6 


MIPS8251 


1.94 


20.33 


1.34 


bulge 


Pre-merger (III) 


2 


MIPS8253 


0.95 


18.64 


3.41 


bulge 


Advanced merger (IV) 


3 


MIPS8308 


0.37 


17.58 


2.80 


disk 


Advanced merger (IV) 


3 


MIPS8311 


1.17 


18.93 


2.99 


disk 


Face-on spiral 




MIPS8315 




19.86 


~ 1.52 b 


ambiguous 


Old merger (V) 


5 


MIPS8325 


0.61 


18.53 


2.27 


bulge 


Triplet 


1 


MIPS8327 


2.44 


20.95 


1.11 


ambiguous 


First contact (II) 


4 


MIPS8328 


1.02 


18.93 


3.10 


disk 


Advanced merger (IV) 


3 


MIPS8342 


1.56 


20.26 


1.21 


bulge 


Advanced merger (IV) 


4 


MIPS8360 


1.50 


18.91 


3.33 


disk 


Advanced merger (IV) 


3 


MIPS8375 


0.87 


18.89 


2.12 


disk 


Pre-merger (III) 


1 


MIPS8387 


0.91 


17.57 


9.11 


disk 


Close pair (I) 


2 


MIPS8388 


1.14 


18.95 


3.38 


bulge 


Advanced merger (IV) 


3 


MIPS8392 


1.90 


21.79 


2.40 


disk 


Faint & Compact 


6 


MIPS8400 


1.51 


18.44 


1.46 


bulge 


Old merger (V) 


5 


MIPS8405 


1.16 


18.66 


1.24 


bulge 


Triplet 


3 


MIPS8407 


0.89 


20.93 


1.32 


ambiguous 


Regular bulge 


6 


MIPS8413 


2.18 


21.57 


1.75 


disk 


Pre-merger (III) 


3 


MIPS8430 


0.67 


18.93 


2.15 


disk 


Close pair (I) 


2 


MIPS8450 


1.00 


19.53 


6.16 


disk 


Face-on spiral 




MIPS8462 


1.01 


18.82 


3.16 


disk 


Close pair (I) 


2 


MIPS8465 




19.33 


~ 0.66 b 


bulge 


Old merger (V) 


4 


MIPS8477 


1.84 


19.44 


1.55 


bulge 


Advanced merger (IV) 


3 


MIPS8479 


1.30 


20.51 


3.25 


disk 


Faint & Compact 


5 


MIPS8493 


1.80 


20.53 


4.16 


disk 


Triplet 


3 
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MIPS ID 


Redshift 


m H 


r i/2 (kpc) 


Dominant 


Morphological 


Confidence 








component 


class 


level 


MIPS8499 


0.60 


17.27 


5.92 


bulge 


Advanced merger (IV) 


2 


MIPS8507 


0.76 


17.98 


4.40 


disk 


Face-on spiral 




MIPS8521 


1.19 


19.58 


3.80 


disk 


Pre-merger (III) 


3 


MIPS8532 


0.86 


18.60 


3.31 


disk 


Face-on spiral 




MIPS8543 


0.65 


18.09 


4.05 


disk 


Face-on spiral 




MIPS15690 


0.85 


18.51 


3.80 


disk 


Edge-on spiral 




MIPS15755 


0.74 


17.66 


4.12 


disk 


Face-on spiral 




MIPS15771 


2.20 


21.07 


1.76 


bulge 


Pre-merger (III) 


4 


MIPS15776 


1.12 


18.65 


2.25 


bulge 


Advanced merger (IV) 


2 


MIPS15840 


2.30 


21.58 


2.20 


disk 


Pre-merger (III) 


4 


MIPS15871 




20.54 


~ 0.64 b 


ambiguous 


Old merger (V) 


5 


MIPS15880 


1.64 


20.84 


3.22 


disk 


Pre-merger (III) 


3 


MIPS15897 


1.62 


18.67 


0.67 


bulge 


Old merger (V) 


5 


MIPS15928 


1.50 


19.59 


2.51 


disk 


Pre-merger (III) 


3 


MIPS15941 


1.23 


20.11 


1.85 


disk 


Faint & Compact 


4 


MIPS15949 


2.12 


20.76 


2.83 


disk 


Advanced merger (IV) 


3 


MIPS15958 


1.97 


21.14 


0.81 


ambiguous 


Old merger (V) 


4 


MIPS15967 


1.30 


18.12 


2.62 


bulge 


Advanced merger (IV) 


4 


MIPS15977 


1.85 


19.71 


2.15 


bulge 


Advanced merger (IV) 


3 


MIPS16037 


1.61 


20.43 


1.86 


bulge 


Regular bulge 


6 


MIPS16047 


0.52 


17.99 


3.46 


disk 


Edge-on spiral 




MIPS16059 


2.33 


20.29 


2.51 


disk 


Pre-merger (III) 


1 


MIPS16080 


2.01 


20.21 


2.36 


disk 


Advanced merger (IV) 


2 


MIPS16095 


1.81 


20.01 


1.75 


disk 


Advanced merger (IV) 


3 


MIPS16099 


0.95 


19.70 


2.02 


bulge 


Advanced merger (IV) 


3 


MIPS16113 


1.90 


20.96 


2.54 


disk 


Pre-merger (III) 


3 


MIPS16118 


2.61 


21.93 


3.50 


disk 


Pre-merger (III) 


4 


MIPS16122 


1.97 


21.40 


2.59 


disk 


Advanced merger (IV) 


3 


MIPS16132 




18.97 


~ 2.22 b 


disk 


Pre-merger (III) 


4 


MIPS16135 


0.62 


19.41 


2.13 


disk 


Edge-on spiral 




MIPS16144 


2.13 


20.36 


2.97 


disk 


Advanced merger (IV) 


3 


MIPS16156 


0.72 


18.16 


3.86 


disk 


Pre-merger (III) 


2 


MIPS16170 


0.32 


17.72 


2.35 


disk 


Edge-on spiral 




MIPS16202 




20.11 


~ 2.11 b 


disk 


Advanced merger (IV) 


3 


MIPS16219 


2.72 


19.68 


1.17 


bulge 


Old merger (V) 


5 


MIPS16249 


0.53 


18.44 


1.80 


disk 


Pre-merger (III) 


2 


MIPS16267 


1.31 


21.11 


1.60 


ambiguous 


Old merger (V) 


4 


MIPS22196 


0.80 


18.73 


0.96 


bulge 


Advanced merger (IV) 


2 


MIPS22204 


1.97 


20.22 


1.05 


bulge 


Pre-merger (III) 


2 


MIPS22248 


1.20 


18.85 


2.54 


disk 


Faint & Compact 


4 


MIPS22277 


1.77 


19.68 


2.02 


ambiguous 


Pre-merger (III) 


4 


MIPS22303 


2.34 


22.60 


1.65 


disk 


Pre-merger (III) 


3 


MIPS22307 


0.70 


18.14 


2.05 


bulge 


Advanced merger (IV) 


2 


MIPS22352 


0.66 


17.67 


3.23 


bulge 


Advanced merger (IV) 


3 


MIPS22371 


1.67 


19.64 


3.07 


disk 


Pre-merger (III) 


3 


MIPS22379 


0.65 


17.85 


3.39 


bulge 


Old merger (V) 


4 


MIPS22397 


1.83 


19.76 


2.78 


disk 


Pre-merger (III) 


3 


MIPS22432 


1.59 


19.80 


4.18 


disk 


Pre-merger (III) 


1 


MIPS22516 


1.35 


19.67 


1.20 


bulge 


Advanced merger (IV) 


3 


MIPS22527 




19.73 


~ 0.84 b 


bulge 


Pre-merger (III) 


4 


MIPS22530 


1.95 


20.70 


3.54 


disk 


Pre-merger (III) 


2 


MIPS22536 


1.59 


19.04 


2.73 


disk 


Triplet 


2 


MIPS22549 


1.05 


20.12 


2.30 


disk 


Faint & Compact 


5 


MIPS22555 


1.88 


20.82 


1.99 


disk 


Pre-merger (III) 


2 


MIPS22557 


0.79 


18.03 


3.17 


bulge 


Advanced merger (IV) 


2 


MIPS22558 


3.20 


21.61 


0.60 


bulge 


Advanced merger (IV) 


3 


MIPS22621 




23.71 


~ 0.67 b 


ambiguous 


Faint & Compact 


6 


MIPS22635 


0.80 


19.36 


1.81 


disk 


First contact (II) 


3 


MIPS22638 


0.98 


19.21 


4.26 


disk 


Pre-merger (III) 


4 


MIPS22651 


1.73 


19.33 


2.50 


bulge 


Pre-merger (III) 


3 


MIPS22661 


1.75 


20.35 


1.70 


bulge 


Advanced merger (IV) 


4 


MIPS22690 


2.07 


21.01 


1.01 


disk 


Pre-merger (III) 


4 


MIPS22699 


2.59 


22.26 


0.85 


ambiguous 


Old merger (V) 


5 


MIPS22725 


1.40 


19.37 


1.57 


bulge 


Old merger (V) 


5 



Note. — The profile of some of our objects is listed as ambiguous. These are objects whose profile is not well constrained either because they 
:e very faint, barely resolved or possess a strong PSF component that hides the underlying galaxy. 

MIPS42 is very close to a diffraction spike in our NICMOS image which causes all of our fits to fail. Nevertheless, it is visually clear that it is 
disk-dominated, "faint & compact" object. 
Assuming 1" = 8 kpc 



